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The  Heuristic  Model  and  Its  Application 
by 

T.  H.  Gawain,  D.Sc. 
and 

J.  K.  Pritchet 
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1.1  INTRODUCTION 


The  Naval  Radiological  Defense  Laboratory,  prior  to  its  disestablish¬ 
ment  in  June  1969,  had  for  some  time  been  pursuing  studies  in  numerical 
hydrodynamics  in  the  field  of  underwater  explosions.  It  is  veil  known 
that  fluid  turbulence  forms  an  important  aspect  of  this  problem.  In  fact, 
the  lack  of  an  adequate  theory  of  turbulence  has  been  an  Important  limita¬ 
tion  on  the  development  of  numerical  hydrodynamics  for  applications  not 
only  to  the  explosion  problem,  but  also  to  fluid  mechanics  problems  of 
every  type.  There  continues  to  be  an  urgent  need  to  overcome  this  obstacle. 

Consequently,  an  agreement  was  reached  between  NRDL  and  the  Naval 
Postgraduate  School  some  time  ago  for  a  joint  attack  on  the  problems  of 
fluid  turbulence.  Two  main  approaches  were  investigated.  The  first  approach 
emphasized  fundamental  understanding  of  the  phenomena  without  undue  concern 
about  developing  computational  methods  for  application  to  practical  problems. 
The  second  approach  emphasizes  the  development  of  a  feasible,  heuristic 
computational  method  at  the  expense,  if  necessary,  of  analytical  rigor. 

This  report  documents  the  excellent  results  that  have  been  achieved  in 
connection  with  this  latter  goal.  Reference  (1),  although  now  largely 
superseded  by  more  recent  work,  documents  some  earlier  effort  along  the 
other  avenue  of  attack.  Further  work  is  continuing  at  the  Postgraduate 
School. 

Unfortunately,  the  Naval  Radiological  Laboratory  was  disestablished 
in  the  midst  of  this  research.  The  abrupt  cessation  of  NRDL  activities 
mcde  it  virtually  impossible  to  complete  certain  details  of  the 
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mathematical  model  In  a  satisfactory  manner.  It  was  therefore  decided  to 
hold  up  distribution  of  a  final  report  until  such  time  as  these  additional 
points  could  be  worked  out.  Meanwhile,  the  press  of  other  sponsored  research 
compelled  the  investigators  to  lay  thia  task  aside  for  a  considerable  time. 

Fortunately,  it  war  eventually  possible  to  complete  and  revise  the 
detalla  mentioned  above.  Consequently,  the  present  report,  which  describes 
the  final  model,  la  now  being  released  for  publication  at  last.  Although 
the  major  part  of  the  work  described  was  accomplished  some  time  ago,  the 
modal  Itself  should  continue  to  be  of  scientific  Interest  for  some  time, 
owing  to  the  basic  and  long  term  nature  of  the  problem  with  which  It  deals. 

As  an  Intereating  hlatorlc  note,  it  seems  worthwhile  to  mention  that 
the  heurlatlc  model  developed  in  this  report  resembles  in  some  respects 
the  approach  first  suggested  by  Prandtl  In  1945,  Ref.  (2).  It  also  has 
some  similarities  to,  and  some  differences  from,  the  more  recent  work  of 
Harlow  and  his  associates.  References  (3)  and  (4) . 

In  Incompressible  flow,  whether  laminar  or  turbulent,  the  equations 
of  momentum  and  continuity,  along  with  the  boundary  conditions,  theoretically 
suffice  to  establish  completely  the  exact  fluid  motion.  If  the  flow  happens 
to  be  turbulent,  however,  the  actual  detailed  motion  becomes  so  complex 
that,  although  it  is  in  principle  determinate,  its  actual  calculation 
would  impose  an  overwhelming  computational  burden. 

The  usual  way  around  this  difficulty  is  to  average  the  equations  of 
momentum  and  continuity.  This  results,  of  course,  in  an  enormous 
simplification.  Unfortunately,  it  also  involves  a  significant  and  irre- 
trlevalbe  loss  of  essential  information.  Consequently,  owing  to  the 
presence  of  the  unknown  Reynolds  stresses  which  are  created  by  this  aver¬ 
aging  process,  the  averaged  equations  of  momentum  and  continuity  do  not 
in  themselves  comprise  a  determinate  set. 
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In  order  to  define  a  determinate  solution,  additional  relations  are 
required  to  fix  the  unknown  Reynolds  stresses.  Unfortunately,  these 
supplementary  relations  cannot  be  established  from  the  original  equations 
by  any  purely  deductive  process.  For  this  purpose,  supplementary  empirical 
hypotheses  are  an  unavoidable  necessity. 

From  another  viewpoint,  it  may  be  stated  that  the  averaged  equations 
of  motion  show  the  effect  of  the  Reynolds  stresses  upon  the  mean  flow. 
However,  the  reciprocal  effect  of  the  mean  flow  upon,  the  Reynolds  stresses 
is  lost  in  the  averaging  process.  Hence  some  adequate  hypothesis  must  be 
found  for  representing  this  relation,  at  least  approximately. 

Fcr  this  purpose,  a  heuristic  approach  which  seems  plausible  is  to 
postulate  a  relation  between  the  Reynolds  stresses  and  the  mean  flow  which 
is  analogous  to  the  relation  that  is  known  to  govern  the  viscous  stresses. 
The  analogue  of  the  ordinary  molecular  kinematic  viscosity  is  the  so- 
called  eddy  kinematic  viscosity.  The  problem  becomes,  therefore,  to 
determine  empirically  the  general  law  which  governs  this  mean  effective 
eddy  viscosity  at  every  space/time  point  in  the  flow  field. 

The  eddy  viscosity  presumably  depends  on  a  number  of  variables,  one 
of  the  most  Important  of  which  is  the  local  kinetic  energy  of  turbulence. 
Therefore,  it  becomes  necessary  to  find  the  space/ time  distribution  of  the 
turbulent  energy.  Fortunately,  the  governing  energy  equation  can  be 
deduced  rigorously  from  the  original  equations  of  motion.  However,  the 
energy  equation  itself  introduces  two  additional  unknowns  which  can  only 
be  approximated  in  the  same  heuristic  and  empirical  fashion  as  the  eddy 
viscosity  Itself.  The  additioral  unknowns  are  the  rate  of  dissipation  of 
turbulent  energy  into  heat,  and  the  rate  of  turbulent  diffusion  of  energy. 

Theory  and  experiment  both  show  that  the  eddy  viscosity,  and  the 
dissipation  and  diffusion  functions  as  well,  depend  not  only  on  the 
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turbulent  energy  itself ,  but  also  on  a  local  length  scale  parameter  which 
can  be  associated  with  each  space/tlme  point  in  the  flow  field.  Von 
Karman  was  perhaps  the  first  to  point  out  how  a  physically  meaningful 
characteristic  length  can  be  defined  in  terms  of  local  space  derivatives 
of  the  mean  velocity  at  any  point  in  the  flow.  In  the  present  paper,  the 
original  approach  of  von  Karman  is  further  developed  and  refined.  It  no* 
takes  into  account  not  only  the  velocity  derivatives  at  the  designated 
point  itself,  but  also  the  values  in  the  general  vicinity  of  the  point. 

By  employing  dimensional  analysis,  and  by  applying  the  available 
experimental  data,  we  finally  obtain  three  empirical  expressions  which 
determine  to  a  reasonable  approximation  the  eddy  viscosity,  the  heat 
dissipation  and  the  turbulent  diffusion,  respectively.  These  expressions 
also  involve  the  turbulent  energy,  the  local  length  parameter,  and  the 
distance  to  the  nearest  fixed  wall  (if  any).  Of  course,  these  empirical 
expressions  are  amenable  to  further  investigation  and  development. 

In  this  way  a  single  consistent  and  determinate  set  of  equations  is 
established  which  applies  in  principle  to  any  incompressible  turbulent 
flow  field.  Only  the  boundary  conditions  differ  for  each  specific 
application. 

The  application  of  this  heuristic  model  is  shown  for  two  widely 
different  cases,  namely,  a  two-dimensional  duct,  and  an  axi-symmetric  jet. 
The  distributions  of  energy  and  Reynolds  stress  are  computed  for  each 
case  and  compared  with  experimental  measurements.  The  agreement  is 
satisfactory. 

These  results,  although  limited  in  extent,  strongly  suggest  that  the 
proposed  heuristic  model  of  fluid  turbulence  is  basically  correct,  at 
least  in  its  main  essentials. 
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This  report  Is  divided  Into  two  parts.  Patt  I  summarizes  the  final 
heuristic  model  and  Illustrates  Its  application  to  two  typical  cases, 
namely,  a  channel  flow  and  a  free  jat.  Part  II  plves  the  corresponding 
detailed  background,  theory  and  rationale.  It  tells  vhy  the  final  model 
Is  Just  what  It  Is.  It  shows  how  the  various  parameters  of  the  model  are 
related  to  fundamental  aspects  of  turbulent  correlation  and  spectra.  An 
original  analysis  is  presented  of  :l»e  general  tensor  relation  between 
stress  and  strain  rate.  This  provides  the  basl3  for  Judging  both  the 
usefulness  and  the  limitations  of  the  eddy  viscosity  concept  as  employed 
In  the  heuristic  model.  The  energy  equation  Is  derived  from  first 
principles,  then  transformed  so  as  to  be  practically  useful.  This  part  of 
the  report  also  raises  various  questions  of  a  fundamental  nature  which 
could  be  fruitful  for  future  research. 
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1.2  SYMBOLS 


a 


A 

b 


A1A2A3^ 

B1B2B3j 


Ref.  Eqs . 

Local  dimensionless  length  scale  in  jet  1.6-31 

Overall  dimensionless  length  scale  of  1.6-28 

cross-section  in  jet 

Cross-sectional  area 

Channel  half  width, or  jet  reference  length  1.6-1 


Functions  in  momentum  equations  for  jet  1.6-9  thru  -14 


2  7  2 

Exponents  in  definitions  of  I  ,  J  ,  X  ,  2.6-14  thru  -19 

and  t 


C  ,  C  Arbitrary  cor  /.ants  in  definitions  of  2. 6-8, -9 

'  X2  and  T2 


D 


Determinant.  Also  diameter  of  jet  2.6-16 


E 


Mean  kinetic  energy  of  turbulence  per  1.3-16,  2.3-5 

unit  mass 


E’ 


Instantaneous  kinetic  energy  of  turbulence  A-12 
per  unit  mass 


H 


Rate  of  dissipation  of  turbulent  energy  2.8-8 
per  unit  mass  into  heat. 


f  Gaussian  function  1.6-6 

f^  Ratio  of  local  length  scale  of  mean  flow  2.5-1 

to  correlation  length 


f 

T 


Ratio  of  local  time  scale  of  mean  flow  2.5-2 

to  correlation  time 


Non-isotroplc  tensor  in  generalized  stress/  2.7-16 
strain  rate  lav 


Denotes  values  of  f..  for  principal  axes 
of  strain  rate  J 
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F 

Dimensionless  stream  function  of  jet 

1.6-1 

F' 

Denotes  derivative 

F,  G,  H. 

Functions  In  energy  equation  for  two- 
dimensional  channel 

1.5-25,  -26-27 

8* 

Parameter  In  stress/strain  rate  lav 

2.7-28 

GV  G2*  G3 

Functions  In  analysis  of  Jet 

1.6-25,  -26,  -27 

H 

Integral  of  weighting  function 

1.3-11 

„(S) 

Integral  of  spatial  weighting  function 

2.5-9 

Bm 

Integral  of  temporal  weighting  function 

2.5-10 

i,  j,  k,  1 

Tensor  Indices 

t2.  j2 

Characteristic  Integrals  for  a  point 

1.3-13  and  -14 

t 

Vcn  Karmen's  mixing  length 

1.5-7,  -8 

(n) 

Superscript  Indicating  realization  number 

N* 

Reynolds  number 

P 

Kinematic  pressure  (pressure  divided  by 
density)  of  turbulent  fluctuations 

• 

P 

Kinematic  pressure  (pressure  divided  by 
density)  of  mean  flow 

q 

Root  mean  square  turbulent  velocity 
fluctuation 

1.3-16 

r,  z,  0 

Cylindrical  coordinates:  radial,  axial, 
and  angular,  respectively 

Rij 

Correlation  tensor 

2.3-6 

R(S) 

Rij 

Spatial  correlation  tensor 

2.3-9 

.(T) 

Rij 

Temporal  correlation  tensor 

2.3-9 

s 

Exponential  constant  In  Gaussian  velocity 
distribution  for  jet  (s  •  0.102) 

1.6-6 

c 

Time 

8 

At 

T 

t 


ij 


r 1 
ij 


U, 


Ui 


U,  V 


-uv  J 

v* 

dv' 

W 

(S) 


(T) 


.(S) 


(T) 


AX 


Tine  increment  2.3-1 

Semi-interval  for  time  averaging  2.2-4 

Viscous  stress  tensor  for  turbulent  A-16 

fluctuations 

Viscous  stress  tensor  for  mean  flow 

Viscous  stress  tensor  for  mean  flow  plus  A-2 
turbulent  fluctuation 

Turbulent  velocity  fluctuation  2.2-2 

Mean  velocity  2.2-1 

Velocity  of  mean  flow  plus  turbulent 
fluctuation 


Radial  and  axial  components,  respectively. 


of  axi-symmetric  mean  flow 

Reference  velocity  of  turbulent  jet  1.6-1 

Kinematic  Reynolds  stress  tensor  (stress  2.7-1 
divided  by  density) 

friction  velocity  1.5-2 

Element  of  volume 

Weighting  function  1.3-10 

Spatial  weighting  function  2.5-4 

Temporal  weighting  function  2.5.5 

Normalized  weighting  function  1.3-12 

Normalized  spatial  weighting  function  2.5-12 

Normalized  temporal  weighting  function  2.5-13 

Position  vector  of  arbitrary  fixed 
reference  point 

Position  vector  of  variable  point 

Displacement  vector  1.3-9 
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.V 


*i  -1 

t 

*,  y.  T  _ 

Cartesian  coordinates 

y 

Distance  to  nearest  nail,  or  distance 
from  centerline  of  jet  or  center  plane 
of  channel 

z 

Axial  coordinate.  Also  nondlmenslonal 
wall  distance  exponent 

1.5-23 

a 

Dimensionless  eddy  viscosity  coefficient 

1.3-5, -18 

P,  6' 

Dimensionless  heat  dissipation  coefficients 

1.3-19,2.8-10 

Y 

Dimensionless  diffusion  coefficient 

1.3-20,2.8-21 

YiJ 

Strain  rate  tensor  of  turbulent 
fluctuations 

A-17 

F1J 

Strain  rate  tensor  of  mean  flow 

1.3-4 

rii*r$2*r!3 

Strain  rate  tensor  components  for  principal 
axes  of  strain  rate 

6iJ 

The  Kronecker  delta  tensor,  with  value  unity 
for  equal  Indices,  but  zero  otherwise 

- 

e 

Effective  eddy  kinematic  viscosity 

2.7-22 

* 

e 

o 

Eddy  viscosity  In  centerline  of  jet 

1.8-16 

cjk*cljkl 

Tensor  versions  of  eddy  viscosity 

2.7-13,-14 

e" 

"Pressure  diffusion"  coefficient 

2.8-17 

e* 

Kinetic  energy  diffusion  coefficient 

2,8-18 

; 

Radial  coordinate 

B- 11 

n 

Dimensionless  distance  from  center  plane 
of  two-dimensional  channel,  or  from 
centerline  of  jet 

1.5-1 

0,  0' 

Angular  coordinate 

0* 

Strain  rate  parameter  In  stress/straln  rste 
law 

2.7-26 

• 

K 

Constant  In  von  Kerman* s  mixing  length 
theory 

1.5-8 

• 

10 

K,  K± 

Wave  number  vector 

1.5-8 

A 

Local  length  scale  of  mean  flow 

1.3-15 

X1 

Overall  length  scale 

1.3-8 

A* 

Correlation  distance 

2.3-11 

XD 

Dissipation  length 

2.8-8 

V 

Molecular  kinematic  viscosity 

§ 

Separation  distance  In  meridional  plane 

B-6 

I 

Summation  sign 

T 

Local  time  scale  of  mean  flow 

2.6-18 

T* 

Correlation  time 

2.3-12 

TiJ 

Kinematic  Reynolds  stress  tensor  (stress 
divided  by  density) 

2.7-1 

i 

ij 

Distortlonal  kinematic  stress  tensor 

2.7-4 

T1J 

Residual  non-lsotropic  kinematic  stress 
tensor 

2.7-19 

<Tlj>v 

Viscous  distortlonal  stress  tensor 

2.7-10 

t'2 

Invariant  of  the  tensor 

2.7-5 

T1,2 

Invariant  of  the  tensor 

2.7-20 

$ 

ij 

Fourier  transform  of  the  correlation 
tensor 

2.4-2 

<l> 

Stream  function  for  jet.  Also  auxiliary 
weighting  function  for  axl-symmetrlc  flows 

1.6-1 

B-  13 

\p& 

Stress  parameter  in  stress/strain  rate  law 

2.7-  29 

u 

Angular  frequency  or  temporal  wave  number 

fi 

Invariant  of  strain  rate  tensor 
representing  generalized  shearing  rate 

1.3-6 

fi' 

Absolute  magnitude  of  the  gradient  of  Q 

1.3-7, 2. 6-5 

Note: 

In  the  jet  analysis,  subscript  1  on  any  quantity 

1.6-33 

denotes  the  values  at  z  ■  1.  It  Is  therefore 
always  a  function  of  n  only. 
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1.3  SlfMMAFY  OF  BASIC  EQUATIONS  OF  THF  HEURISTIC  MODEL 


The  purpose  of  this  section  Is  to  summarize  the  lasic  equations 
which,  taken  together,  constitute  the  proposed  heuristic  model  of  fluid 
turbulence.  A  few  explanatory  comments  arc  offered  to  identify  these 
various  equations  and  to  show  their  general  relation  tr  the  overall  theory, 
but  no  detailed  justification  is  given  in  this  section;  such  details  are 
Included  in  other  parts  of  this  report.  Also,  the  equations  in  this 
section  are  sthted  in  three  dimensional  cartesian  form.  The  reductions 
to  forms  specifically  applicable  to  plane  or  axl-symmetrlcal  mean  flows 
are  treated  elsewhere. 

An  incompressible  mean  flow  field  must  everywhere  satisfy  the 
continuity  condition 


-  0 


0.3-1) 


and  the  ensemble  averaged  Navler  Stokes  equations  of  motion 


f!!il  +  JLfr  „  1 

f3Ul  1 

.  fw_)  +  fllul 

Ist  1  llJ. 

U«,l  l»*j  1 

as  well  as  the  particular  boundary  conditions  appropriate  to  the  specific 
flow  field  under  consideration. 

However,  because  of  the  occurence  of  the  unknown  Reynolds  stress 
F.qs.  (-1)  and  (-2)  do  not  in  themselves  constitute  a  determinate  set 
of  equations.  In  order  to  define  a  determinate  solution,  additional  relations 
are  required.  Purely  deductive  theory  alone  does  not  provide  these  relations. 
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WsW)W»^niy>«>c< wfww- 


Hence  appropriate  supplementary  hypotheses  must  be  tentatively  postulated, 
subject  to  eventual  verification,  modification  or  disproof  by  experimental 
evidence . 

For  the  reason  discussed  at  length  elsewhere,  we  postulate  that 
the  Reynolds  stresses  can  be  adequately  related  to  the  strain  rates  of  the 
mean  flow  through  a  law  which  has  been  deliberately  somewhat  oversimplified 
to  the  form 


1  2.  .  _ 

Tlj  •  *  3  ’  4iJ  +  cri3 

2 

In  this  equation,  q  is  the  mean  square  velocity  of  the  turbulent 
fluctuations,  c  is  the  so-called  eddy  kinematic  viscosity,  and  is  the 
strain  rate  of  the  mean  flow  as  defined  by  the  expression 


ij 


axjj 


(-4) 


Unfortunately,  the  eddy  viscosity  e  as  defined  above  is  not  a  simple 
property  of  the  fluid  Itself,  but  rather  a  complex  property  of  the  turbulent 
flow  field.  Physical  and  dimensional  reasoning  leads  to  a  supplementary 
postulate,  namely  that  c  is  expressible  in  the  following  form. 


■  a  (-5) 

Here  X  represents  a  local  length  scale,  defined  below,  which  varies 
from  point  to  point  in  the  field,  but  which  is  everywhere  related  in  a 
consistent  way  to  the  local  mean  velocity  distribution.  The  dimensionless 
coefficient  a  is  a  slowly  varying  universal  function  not  predictable  from 
theory,  but  determinable  from  experimental  data.  This  function  is  specified 
later  in  Eq.  (-18) 
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To  establish  the  quantity  A  at  any  point  5,  we  proceed  as  follows. 
We  first  compute  the  following  mean  flow  variables  throughout  the  field, 
namely. 


(-6) 

(-7) 


Then  we  compute  the  reference  quantity 


X2  . 

1  / (flft ' ) 2dv 


(-8) 


o 

where  the  volume  Integrals  extend  over  the  entire  field.  Of  course  A^  Is 
therefore  a  characteristic  constant  for  the  field  as  a  whole  and  Is 
Independent  of  but  It  may  vary  as  a  function  of  time  If  the  mean  flow 
be  unsteady. 

The  calculation  continues  as  follows.  All  Integrals  extend  over 
the  entire  field  of  flow.  Let 

x  -  coordinate  of  the  fixed  point  at  which  A  Is  being  evaluated 
S'  •  coordinate  of  variable  point  at  which  Integrands  below  are 
evaluated 


AS  •  (S'  -  ID  •  separation  variable 


(-9) 


Ax*AS 


weighting  function 


(-10) 


H  ■  /Wdv'  »  normalizing  integral 
W 

w  ■  —  ■  normalized  weighting  function 


(-11) 

(-12) 
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2  A 

I  ■  /wfi  dv'  “  first  characteristic  integral 
2  2 

J  ■  /v(fifl')  dv  ■  second  characteristic  integral 

2  2 

Of  course  the  Integrals  T  and  J  are  functions  of  x  and,  if  the 
mean  flow  is  unsteady,  of  time  t  as  well. 

The  required  length  scale  A  at  any  point  x  is  finally  given  by  the 

ratio 

x2.i! 

A  o 

The  above  calculation  sequence  establishes  not  only  the  desired 
local  length  scale  X  but  also  the  useful  auxiliary  quantity  J.  It  turns 
out  that  J  is  needed  below  in  connection  with  the  energy  equation.  The 
rationale  involved  is  fully  explained  elsewhere  in  this  report. 

If  the  mean  flow  happens  to  be  either  plane  or  axi-symmetric,  all 
of  the  above  volume  Integrals  can  be  reduced  to  corresponding  surface 
integrals.  The  reduction  for  the  axi-symmetric  case  is  rather  interesting; 
it  is  given  in  Appendix  B. 

If  the  flow  field  happens  to  be  self-similar,  the  above  calculations 
need  be  made  on  only  one  representative  cross-section.  Corresponding 
results  for  other  cross-sections  are  then  found  simply  by  applying  the 
proper  scaling  rules.  If  the  flow  field  is  not  self-similar,  all  points 
must  be  computed  Independently.  If  the  flow  field  is  steady,  the  calcula¬ 
tions  can  be  made  once  and  for  all.  If  it  is  unsteady,  they  must,  of 
course,  be  repeated  at  appropriate  intervals  of  time. 

Once  the  distribution  of  A  has  finally  been  established,  it  is  then 
necessary  to  determine  the  distribution  of  the  turbulent  energy  E,  or 
equivalently,  of  the  root  mean  square  turbulent  velocity 
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(-13) 

(-14) 


(-15) 


q  -  v^E 


(-16) 


This  Is  governed  by  the  energy  equation 


Of  course  the  eddy  viscosity  z  in  the  energy  equation  is  the  same 
as  that  earlier  defined  in  Eq.  (-5)  in  connection  with  the  stress/strain 
rate  law.  If  desired,  Eq.  (-5)  may  be  used  to  eliminate  c  from  the 
energy  equation.  In  any  case,  the  coefficient  a  which  governs  c  is  fixed 
by  Eq.  (-18)  below. 

Each  term  in  Eq .  (-17)  has  dimensions  of  time  rate  of  change  of 
energy  per  unit  mass.  The  four  terms  on  the  right  represent,  respectively, 
generation  of  turbulent  energy,  dissipation  to  heat,  convection  by  the 
mean  flow,  .and  turbulent  diffusion.  Molecular  diffusion  is  negligible. 

The  dimensionless  coefficients  8  and  y  which  occur  in  tha  energy 
eouatlon,  like  the  dimensionless  coefficient  a  introduced  earlier,  are 
assumed  to  be  well  behaved  universal  functions.  However,  they  are  not 
determinable  from  theory,  but  must  be  established  empirically. 

On  the  basis  of  the  work  done  to  far,  the  following  expressions 
seem  to  provide  satisfactory  agreement  with  the  available  data. 


(-18) 


(-19) 


(-20) 
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In  these  expressions,  y  Is  the  distance  to  the  nearest  fixed 
boundary.  As  the  boundary  is  approached  y  approaches  zero,  but  it  turns 
out  that  under  these  conditions,  the  length  parameter  k  also  approaches 
zero  simultaneously,  so  that  the  ratio  ^  remains  finite.  In  fact,  it 
happens  that  at  the  wall  Itself,  f^]  equals  unity  (not  zero).  In  any  case. 


it  is  stipulated  in  connection  with  the  above  formulas,  that  %  shall 


arbitrarily  be  assigned  a  lower  limit  of  unity.  On  the  other  hand,  far 
from  any  fixed  boundary,  or  in  the  absence  of  fixed  boundaries,  ve  set 
y  *  •,  whereupon  a,  {?,  and  y  reduce  to  three  simple  constants. 

The  proposed  heuristic  model  of  fluid  turbulence  consists,  therefore, 
of  tne  foregoing  equations  (-1)  through  (-20)  inclusive,  plus  the  particular 
boundary  conditions  appropriate  to  the  specific  problem  under  consideration. 
This  system  of  equations  theoretically  defines  a  determinate  solution  for 
the  general  case  of  Inhomogeneous  and  non-stationary  fluid  turbulence.  It 
fixes  the  mean  flow  field,  the  distribution  of  turbulent  energy  and  the 
Reynolds  stresses. 

The  results  obtained  so  far  tend  to  support  the  proposed  model,  but 
further  work  is  necessary  to  ascertain  more  fully  its  adequacy  and  range  of 
applicability.  Various  features  of  the  model  are  readily  amenable  to 
modification  if  necessary,  so  as  to  improve  its  range  and  accuracy. 

An  important  merit  of  the  model,  if  it  continues  to  prove  successful, 
is  that  it  reduces  a  vide  range  cf  dissimilar  flow  geometries  and  flow 
conditions  to  a  single  unified  theory. 


1.4  CENEPAI  b'ETHOP  OF  SOLUTION 


Although  the  equations  (-1)  through  (-20)  of  the  previous  section, 
along  with  the  appropriate  boundary  conditions,  theoretically  suffice  to 


determine  all  characteristics  of  the  flow  field,  the  actual  solution  of 
this  system  of  equations  for  any  non-trivial  case  is  hy  no  means  simple. 

Let  us  consider  first  the  general  case  of  an  unsteady  mean  flow,  with  the 
initial  velocity  distribution  of  time  t  ■  0  arbitrarily  prescribed. 

Of  course,  the  prescribed  initial  velocity  distribution  must  satisfy 
the  continuity  condition,  Eq.  (-1).  This  offers  no  particular  difficulty. 

In  the  usual  case  in  which  the  mean  flow  is  either  plane  or  axl -symmetric, 
the  continuity  condition  may  be  satisfied  identically  by  expressing  the 
velocity  components  in  terms  of  a  stream  function.  Let  us  suppose  that 
this  has  been  done,  and  that  the  corresponding  mean  flow  velocity  components 
at  time  t*0  have  been  calculated.  The  subsequent  calculation  procedure  may 
now  be  described  in  somewhat  general  terms  a6  follows. 

The  next  step  is  to  perform  the  sequence  of  calculations  summarised 
in  Eqs.  (-4)  through  (-15).  In  this  way  is  obtained  the  local  length 
scale  A  at  every  point  in  the  field,  and  the  quantity  J. 

Kext  we  must  solve  the  energy  equation  (-17)  for  the  mean  square 

turbulent  velocity  q  at  every  field  point.  For  this  purpose  we  substitute 

F.q.  (-5)  into  (-17)  and  use  Eqs.  (-18),  (-19),  and  (-20)  to  evaluate  the 

three  empirical  coefficients  a,  f,  and  v.  (In  some  cases  the  actual 

2 

distribution  of  the  turbulent  energy  E  -  ^  may  be  known  from  experimental 
measurement.  If  so,  we  have  at  this  point  one  important  comparison  between 
theory  and  experiment.) 

Once  estimates  of  A  and  o  are  available  from  the  foregoing  calculations, 
it  becomes  a  simple  matter  to  solve  Eq.  (-5)  for  the  eddy  viscosity  e,  and 
Eq.  (-3)  for  the  Feynolds  stresses. 
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The  F.eynolds  stresses  obtained  in  this  way  may  now  be  substituted 
back  into  the  Navier  Stokes  relations,  Fes.  (-2).  By  taking  the  divergence 
of  Eqs.  (-2),  we  find  that 


(-A.1  .  -  B 

^3x13x  J  3x13x  1  i 


j  + V 


(1.4-1) 


Since  mean  velocity  and  stress  are  at  this  point  known,  the  above  result 
can  he  solved  for  the  pressure  distribution,  subject  to  the  specified 
boundary  conditions  on  pressure.  This  part  of  the  solution  involves  a 
lengthy  iteration  procedure,  but  is  straightforward  in  principle. 

Once  li . ,  t.  ,  and  P  are  all  established,  we  can  revert  to  Eq.  (-2) 

*  .M*  . 


and  solve  for  the  time  rate  of  chanpe  of  velocity  at  every  field 

point.  By  integrating  forward  through  time  over  a  small  interval,  we 
can  find  the  new  mean  velocity  distribution  at  the  later  time. 

The  above  procedure  may  be  repeated  as  often  as  desired,  thus 
tracing  the  gradual  evolution  of  the  flow  field  over  time.  If  the  boundary 
conditions  happen  to  be  stationary  in  time,  the  entire  flow  field  will 
evolve  toward  a  stationary  state,  thus  defining  the  solution  of  the 
corresponding  steady  state  problem. 

The  above  process  of  reaching  a  steady  state  asymptotically  in 
time  can  be  quite  expensive  in  terms  of  computational  effort  and  time 
required.  It  is  often  possible  to  shorten  or  abridge  the  above  calculation 
in  various  ways,  depending  on  the  situation.  For  example,  instead  of 
substituting  the  Reynolds  stresses  predicted  by  Fo.  (-3)  into  the  Navier 
Stokes  law  F.q.  (-2),  we  can  solve  the  latter  for  the  Reynolds  stresses 
'“required"  in  order  to  maintain  the  prescribed  mean  motion  in  a  steady 
state.  These  "required"  stresses  can  then  be  compared  with  the  stresses 
’predicted"  from  Eq.  (-3).  If  the  agreement  between  the  two  sets  of 
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stresses  is  deemed  satisfactory,  this  can  be  taken  as  a  verification  of  the 
adequacy  of  the  given  velocity  distribution  as  a  estimate  of  the  steady 
state.  This  is  the  procedure  that  has  been  followed  in  our  sample  problems, 
and  it  has  proved  entirely  satisfactory  for  the  intended  purpose.  In 
this  case  the  purpose  has  been  to  demonstrate  the  degree  of  agreement  of 
our  heuristic  model  with  experimental  observations.  The  agreement  is 
considered  satisfactory  for  both  cases  analyzed. 


1.5  APPLICATION  TO  A  TWO-DIMENSIONAL  CHANNEL 

Consider  the  steady,  incompressible,  turbulent  flow  through  a  two- 
dimensional  channel  of  uniform  h  .sight  2b  as  shown  in  Fig.  1.5-1.  Let  y 
be  distance  from  the  channel  center-line.  Choosing  b  as  the  reference  unit 
of  length,  the  corresponding  dimensionless  distance  becomes 

J  -  n  (1.5-1) 

Let  the  kinematic  shear  stress  at  the  wall  be 

-  v*2  (-2) 

where  v*  is  the  so-called  friction  velocity.  Without  loss  of  generality, 
we  can  always  choose  reference  units  of  time  such  that 

■  v*2  ■  1  (-3) 

The  kinematic  shear  stress  i  varies  linearly  across  the  channel  so 
that  at  any  location  n  we  have 

t  -  n  (-4) 

According  to  our  present  notation,  if  l'  is  the  mean  velocity  at  any 

n,  then 
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-  U'  -  fi 


(-5) 


dU 


dn 


{d2V 


dn 


-  u':  -  £>• 


According  to  the  von  Karman  mixing  length  theory 
2  2 

t  -  l  n 

where  the  mixing  length  l  is  given,  in  our  notation,  by 

.2 


.  2  2  n 

l  ■  K  — 

fi 


,2 


(-6) 


(-7) 


(-8) 


and  where  k  is  a  constant  which  remains  to  be  determined. 

Upon  combining  Eqs.  (-4),  (-7),  and  (-8),  taking  the  square  root 
of  both  sides  ard  rearranging,  we  obtain  the  basic  relation 


P'  K 
IT  /iT 

This  equation  is  readily  integrable  and  gives  the  well  known  von 
Karman  velocity  profile  which  in  the  present  notation  becomes 


(-9) 


-f  -  U  -  U0  +  ^  Un(l  -  /n)  +  /n)  (-10) 

It  is  known  that  this  results  gives  satisfactory  agreement  with  experiment 
if  we  place 

<  -  0.36  (-11) 

The  resulting  quantities  of  interest  in  the  present  analysis  are 


fi  -  l”  »  - 


2k(1  -  /n) 


and 


fi'  -  U"  ■  - 


Ate  n(l  -  /n)' 


(-12) 


(-13) 
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We  next  require  the  "averaging  distance"  Xj  defined  hy  the  equation 


C  -  Urn 

,  r  i 

V1  Jo 


Pi  4 

J  ft  dn 


(-14) 


(ftft')2dn 


Upon  substituting  expressions  (-12)  and  (-13)  into  (-14),  we  find 
that  In  the  limit  as  rij  approaches  unity,  both  Integrals  approach  Infinity, 
hut  the  above  ratio  of  these  Integrals  approaches  zero.  Hence  we  obtain 
the  unusual  result  that,  In  this  case, 


X2  -  0 

Now  the  local  length  scale  X  may  be  computed  from  the  relation 


(-15) 


1  AX’ AX 


e  X 


1  J**  . 

fl  dn 


Mm  ^ 

^•*0  f1-  p' 

I  ®  2 

Jo  1  (nn')  dn’ 


This  result  Is  remarkable!  It  shows  that  our  characteristic 
length  X  is  In  this  case  proportional  to  the  classical  von  Karman  mixing 
length  4.  Form  Eqs.  (-8),  (-12),  (-13)  and  (-16)  we  find  that 


(-16) 


X  ■  ~  *  ~r  •  2/n(l  -  /n) 

Similarly,  because  of  the  exceptional  result  (-15) ,  the  Integral 

2 

for  J  reduces  simply  to 

( 1  A  x.  A  x 


e  1  (ftft')2dn' 


J  -  Mm 

xx-*o 


1  AX*  AX 
'  >1 

e  1  dn' 


-  (ftft')2 


(-17) 


(-18) 


o 


Therefore  again  using  (-12)  and  (-13)  we  find  that 


J  ■  fin'  -  — - — l - 

8</n(l  -  /n)3 


(-19) 


The  three  empirical  functions  a,  B,  and  y  may  now  be  expressed  in 


the  form 


■  a  -  0.065(1  +  e"z  } 


(-20) 


j  “  3.7(1  +  e"Z  } 


(-21) 


Y  -  1.4  -  0.4  e~z 

where 

M/H2 

With  the  functions  fi.  A,  J,  a,  g,  and  Y  now  knownf  wp  can  proceed 
to  solve  the  energy  equation.  The  unsteady  term  and  the  convective  term 
both  vanish  from  this  equation  in  the  present  instance  and  we  obtain 


(-22) 


(-23) 


oXVSEfi2  -  B(2E)7/£J1/3  +  jL  «YX/a[gjl-  0 


(-24) 


Energy  E  is  now  the  only  unknown.  To  simplify  the  equation  further. 


we  introduce  the  following  notation. 

aAfi2  -  F(n) 

8J1/3  -  G(n) 
ay\  ■  H(n) 

Of  course,  these  three  functions  are  now  known.  Therefore  the  energy 
equation  reduces  to 


(-25) 


(-26) 

(-27) 
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G(2E)7/6  +  -  0 


(-28) 


This  equation  cannot  be  solved  analytically,  but  it  can  easily  be 
expanded  and  rearranged  for  numerical  solution  by  finite  differences. 

The  desired  boundary  conditions  are 


At  n 

V' 


0  i¥\  - 


E'  -  0 
o 


(00 


PC  At  n  ■  l  F.^  -  0 


(-29) 


(-30) 


This  last  OOridltion  cannot  be  satisfied  directly.  However,  vA  can  integrate 

from  the  centerline  out,  using  a  trial  value  of  the  centerline  energy  E  . 

° 

By  trial  and  error,  a  value  of  E  can  then  be  found  which  will  in  fact 

°  < 

I 

yield  the  required  condition  £.  ■  0  at  the  wall. 

1  V. 

This  is  the  method  which  was  used  to % compute  the  energy  distribution 
shown  by  the  solid  line,  in  Fig.  1.5-2.  It  may  be  seen  that  the  curve 

obtained  in  this  way  is  in  reasonably  good  agreement  with  the  experimental. . 

> 

data  of  Relchart  and  of  Laufer,  Refs.  5  and  6. 

Once  the  energy  distribution  is  known,  the  dimensionless  Reynolds 

A  v' 

p.r 

stresses  may  be  computed  from  the  relation 


-uv  «  t  ■ 


(-31) 


Except  very  near  the  walls,  the  above  Reynolds  stress  should  agt'ee 
with  the  total  shear  stress  defined  by  Eq.  (-4). 

If  we  wish  to  switch  back  to  absolute  units,  Eqs.  (-4)  and  (-31) 
reduce  respectively  to 


T  V 

* 


Total  Stress 


(-32) 


uv  fx)  fyjf]  ft 


Reynolds  Stress 


(-33) 
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FIGURE  1.5-2  TURBULENT  ENERGY  DISTRIBUTION  IN  A 
TWO  DIMENSIONAL  CHANNEL 


% 
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In  Fig.  1.5-3,  the  solid  line  shows  the  above  heuristic  estimate  of 
the  Reynolds  stress,  computed  for  a  Reynolds  number  of  10,000.  The  dashed 
straight  line  represents  the  true  total  stress,  which  should  agree  closely 
with  the  Reynolds  stresr  everywhere  except  in  the  immediate  vicinity  of 
the  wall.  The  agreement  is  seen  to  be  good;  the  computed  results  are  on 
the  whole  more  accurate  than  are  the  experimental  measurements  from  Reference 
(6),  as  shown  on  this  figure. 

We  conclude  that  this  example  tends  to  substantiate  the  proposed 
heuristic  model  of  turbulence. 

1.6  APPLICATION  TO  AN  AXI-SYMMETRIC  JET 


Consider  a  free  turbulent  jet  discharging  into  a  quiescent  atmosphere 
as  shown  in  Fig.  1.6-1.  The  radial  and  axial  coordinates  are  r  and  s  and 
the  corresponding  velocity  components  are  U  and  V,  respectively. 

Experimental  results  indicate  that  the  velocity  profiles  of  the 
mean  flow  at  various  cross-sections  are  self-similar.  Specifically,  if 
z  be  measured  from  a  suitable  virtual  origin,  the  stream  function  ip  can 
be  reduced  to  the  form 


ip  •  U  bzF(n) 
o 

where 

r 

n  -- 

and  where  U  and  b  are  constants  of  the  jet. 
o 

Consequently  the  velocity  components  of  the  mean  flow  become 


(1.6-1) 


(-2) 


(-3) 
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FIGURE  1.6-1  THE  AXI- SYMMETRIC  JET 


(-4) 


where 


F' 


'3F' 


(“5) 


The  known  experimental  results  also  indicate  that  the  generalized 
velocity  distribution  through  the  jet  Is  very  well  approximated  by  the 
simple  expression 

2 

_  n_ 

F'  .  s2 

“  -  f  -  e  (-6) 

Where  s  la  a  characteristic  constant  for  a  turbulent  jet*  It  has  the  known 
value  8  -  0.102. 


Furthermore,  It  Is  useful  and  convenient  to  adopt  units  of  length 

* 

and  time  such  that,  when  expressed  in  these  units,  IJo  •  1  and  fc  ■  1. 

Of  course,  there  Is  no  loss  In  generality  involved  In  this.  It  amounts 
merely  to  non-dlmenslonallzlng  all  quantities  with  respect  to  Uq  and  b  as 
reference  parameters.  Hence  symbols  like  e,  p,  u,  v  and  so  on  will  now 
represent  the  corresponding  dimensionless  versions  of  kinematic  viscosity, 
kinematic  pressure,  velocity  components  and  the  like. 

The  next  step  In  the  analysis  is  to  combine  the  momentum  equations 
(-2)  and  the  stress  law  (-3)  of  section  1.3.  The  results  may  he  transformed 
to  the  z,  n  coordinate  system  defined  above.  The  foregoing  dimensionless 
functions  c  and  p  may  be  Introduced  therein,  and  the  results  reduced. 

This  reduction,  although  very  lengthy,  is  straightforward.  In  this  way 
the  following  results  are  obtained. 

np'  +  2p+  +  A2(n)e  ■  A3(n)  (-7) 

-p'  +  0  +  B1(n)e'  +  B2^n^c  "  ^(h)  (-8) 
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The  A's  and  B's  are  known  functions  defined  as  follows,  where  f  is  the 


Gaussian  function  defined  by  Eq.  (-6). 

V> '  (i  -  h 

v  s  s 


:<">  ■  [(2  -  ?)  -  (s  -  7]  ¥  +  ¥]  * 


(-10) 


/3(n)  -  f (1  -  2f) 


‘  ('  '  72  '  +  (1  ' 


»2(n,  -  f  (l  -  $  -  (3  -  ^ 


/  2  /  M 

*n  .  An_  f 

T+jT  nf 

8  8  J 


(-11) 


(-1?) 


(-13) 


o  S2'  o 

B3(n)  -  -nf  +  n  +  f(l  -  f)  -  — 3  (1  -  f)z  (-14) 

'  '  4n 

Theoretically,  it  is  possible  to  eliminate  p  and  p'  between  the  two 
momentum  equations  (-7)  and  (-8).  This  yields  a  single  vorticity  eauatlon 
of  second  order  in  e  which  can  be  solved  for  the  unknown  eddy  viscosity 
distribution  e(n). 

However,  it  is  known  from  experiment  that  the  terms  in  p  and  p'  are 
negligible  in  Eqs.  (-7)  and  (-8).  Moreover,  the  terms  in  Eq.  (-8)  are 
negligible  in  comparison  with  those  in  Ec.  (-7).  We  can  take  advantage  of 
this  information  to  simplify  the  solution  for  e(n).  We  simply  discard 
Eq.  (-8),  and  drop  the  terms  in  p  and  p'  from  Lq.  (-7).  The  resulting 
equation  is  now  only  of  first  order  in  e.  It  may  be  written  in  the  form 


*_]  .  U  -  2-1  +  *> 

2  2  2  * 
s  s  ;  s  _ _s 

f.  2_  2n21 

I  s*  s2Jn 


:i  -  2f) 


fi  1_  lull 
l1  "  ?  “  2  n 

^  fi  fi 


(-15) 


Solution  of  this  equation  analytically  Is  not  feasible,  but  solution 
by  numerical  integration  is  simple. 

bote  first  that  at  the  centerline  n  ■  0,  this  eauatlon  gives  at  once 


TlT 

L2 


-  0.00260 


(-If) 


This  result  provides  the  boundary  value  required  for  the  numerical  integration. 

It  1 8  also  of  interest  to  note  that  for  very  large  values  of  n, 
tq.  (-15)  simplifies  to  the  following  approximation,  namely, 


(-17) 


The  solution  is 


E 


(-18) 


The  true  numerical  solution  approaches  the  above  analytical  curve  asymptoti¬ 
cally  at  large  values  of  r,. 

With  the  dimensionless  eddy  viscosity  distribution  c(n)  Vnown,  it 

becomes  a  simple  matter  to  find  the  corresponding  dimensionless  Reynolds 

shear  stress  t  from  the  relation 
rz 


rrz  ■  —  *  ‘[(£)  +  (£)] 

By  making  use  of  Eqs.  (-3)  through  (-6)  this  result  can  be  reduced  to  the 
following  form  for  computation 


The  shear  stress  computed  from  this  enuation  is  shown  by  the  dashed  line 
in  Fig  (1.6-2) . 
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To  form  the  mean  flow  parameters  required  for  the  heuristic  model, 
we  note  that,  in  cylindrical  coordinates, 


2 


+  2 


'3  V'  2 


?z 


+  2 


I2 + M +  (£)] 


and 


The  change  to  the  new  coordinates  n,  z  is  quite  lengthy,  but  not 
difficult  otherwise.  The  results  finally  obtained  can  be  summarized  in 
the  following  form. 


fl2  -  \  G^n) 
z 


(ra 


,)2  "  ~io  [c,2(r,)  +  G*(n) 

9 


■] 


2  3 

where  the  G’s  are  known  functions  defined  es  follows. 
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Note  chat  these  three  functions  are  all  finite  on  the  axis  p  ■  0 
and  that  they  all  vanish  at  infinity. 

To  find  the  averaging  distance  X^  in  a  non-dimensional  form  appropriate 
to  this  problem,  we  evaluate  the  expression 


,\^2  r  a4rdr  r  G?(n)ndn 

I  J.|  2  _  *  0  '  0  1  ,  _ 

T-  “  ai  "  - 5 -  '  - 2 - 2 -  (’2 

1  >  /0(«2')  rdr  Jq  (c|(n)  +  C2(n)]ndn 

It  should  be  observed  that  both  integrals  in  this  last  equation  are 
evaluated  from  n  •  0  to  n  »  »  over  the  surface  z  •  constant,  and  that 
in  this  case  the  quantity  z  cancels  identically  from  the  right  sice  of 
the  equation.  There  is  no  integration  with  respect  to  z.  The  dimensionless 
averaging  distance  a^  is  constant  and  Independent  of  z. 

The  Integrals  Involved  in  Eq  (-28)  can  easily  he  evaluated  numerically, 
thus  fixing  a^.  It  is  found  that  a^  is  roughly  the  same  as  s,  the  jet 
width  parameter;  this  seems  reasonable. 

2 

A  similar  procedure  is  involved  in  determining  the  integrals  I  and 

2 

J  and  the  resulting  local  length  scale  parameter.  Also  the  effects  of 
polar  symmetry  must  be  taken  into  account  in  the  manner  of  Appendix  B. 

Again  the  integrations  are  from  n  ■  0  to  n  ■  °>  over  the  surface  z  ■ 
constant,  and  again  the  variable  cancels  from  the  final  result.  Thus 
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Hence  the  dimensionless  local  length  scale  becomes 


M2  2,  ,  1  f I'l 

tJ  “a(n)‘7l7i 


I*(n) 
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2  2 

The  Integrals  I^(n)  and  J^(n)  must  of  course  le  evaluated  numerically. 
It  turns  out  that  the  local  length  scale  a(n)  ohtalred  in  this  way  Is  a 
slowly  varying  function.  The  value  remains  nearly  constant  at  about  0.1 
across  most  of  the  jet,  and  declines  only  slowly  at  relatively  large 
distances  from  the  centerline. 

The  energy  equation  for  the  axl-symmetrlc  turbulent  jet  may  be 
vritten,  first  In  cylindrical  coord!;  -tes. 


[3E 
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(-32) 
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The  various  functions  which  appear  In  the  energy  enuation  can  now 
be  represented  as  functions  of  the  similarity  coordinates  z  and  r,  as 
follows 
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Upor  substituting  the  expressiona(-33)  Into  the  energy  equation,  ve 
obtain , 

^  .  <2Ei>7/6jiin  i  »  .  „  t  ,  a  \\h] 
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(-34) 


This  result  confirms  the  fact  that  all  terms  of  the  enerpy  equation 
>4 

are  proportional  to  z  and  that  the  above  energy  law  Is  Indeed  compatible 
with  the  self-similarity  of  the  overall  solution.  In  order  to  expand 
this  so  as  to  be  able  to  eliminate  z  from  the  result,  we  make  use  of  the 
relation 
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In  this  way  we  obtain  the  general  result 
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Considerable  simplification  can  be  achieved  with  little  loss  of 
accuracy  by  neglecting  the  effect  of  the  slight  divergence  of  the  jet 
on  the  energy  equation.  This  amounts  to  the  assumption  that 


h[  Jr®°  <-”> 

In  this  case  the  energy  equation  simplifies  to 

-  8<2ei>7/V/3  -  i  iWi> +  i  0  <-3S> 

To  complete  the  solution  it  is  necessary  to  knov  the  empirical 
functions  a,  8,  and  y  as  defined  by  Eqs.  (-18),  (-19),  and  (-20)  of 
Section  1.3.  In  this  case,  however,  there  are  no  fixed  walls  Involved, 
and  these  quantities  become  constants,  namely, 
a  -  0.065  , 


Y  ■  1.4 


The  eddy  viscosity  c(n)  which  occurs  in  .he  energy  equation  may  be 
found  in  either  of  two  ways.  );„atly,  we  may  use  the  solution  previously 
found  from  the  momentum  analysis,  by  numerical  integration  of  I:q.  (  -15). 
Secondly,  we  may  use  the  basic  postulate  of  the  heuristic  model,  that  is, 


e  ■  oXTSlE  ■  aaJSE^ 


(-40) 


This  second  alternative  represents  a  somewhat  more  severe  test  of 
the  heuristic  model  and  was  adopted  for  the  numerical  example  here 
discussed.  Thus  substituting  Eq.  (-40)  into  (-38)  gives 


aa^fij  -  E(2E1)7/6J11/3 


7  h[nW  * 


•  0 

(-41) 
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The  energy  function  E^(n)  is  the  only  unknown  remaining  in  this 
equation.  The  boundary  conditions  are 


At 

n  ■  0 

■3Et' 

[an 

-  0 

(-42) 

At 

n-H» 

El~ 

(-43) 

Eq.  (-41)  must  be  solved  numerically.  It  happens  to  be  convenient 
for  this  purpose  to  use  the  substitution 

2E,  -  q*  (-44) 


The  equation  is  then  expanded  and  the  derivatives  replaced  by  appropriate 
finite  difference  approximations  in  a  straightforward  manner. 

The  boundary  condition  (-43)  at  infinity  must  be  satisfied  indirectly. 


Actually,  the  integration  starts  at  the  centerline  using 


M 

W. 


0  and 


a  trial  value  of  qo<  A  value  of  qQ  is  found  by  iteration  which  yields 
q*  0  at  some  large  value  of  n.  In  this  case  the  value  n  ■  0.25  may  be 
regarded  as  "large.1'  It  is  found  that  the  precise  positloi  of  the  outer 
boundary  makes  virtually  no  difference  in  the  results  as  long  as  it  is 
greater  than  about  0.20. 

The  energy  distribution  found  in  this  way  Is  shown  by  the  solid  line 
in  Fig.  1.6-2.  The  agreement  with  the  experimental  results  of  Corrsin, 

Ref.  (7),  Laurence,  Ref.  (8),  is  seen  to  be  satisfactory,  especially  in 
view  of  the  degree  of  scatter  in  the  data  points. 

Once  the  energy  distribution  is  known  the  eddy  viscosity  is  fixed 
by  Eq.  (-40)  and  the  corresponding  shear  stresses  by  Eq.  (-20).  These  are 
shown  by  the  solid  line  in  Fig.  1.6-3.  The  dashed  line  shows  the  corres¬ 
ponding  shear  stresses  as  computed  from  the  momentum  equation.  Data  points 
are  from  Corrsin,  Ref.  (9). 
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FIGURE  1.6-2  TURBULENT  ENERGY  DISTRIBUTION  IN 
A  CIRCULAR  JET 
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FIGURE  1.6-3  REYNOLDS  SHEAR  STRESS  DISTRIBUTION 
IN  A  CIRCULAR  JET 
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If  the  heuristic  model  were  completely  correct,  and  If  the  assumed 
Gaussian  velocity  profile  were  exact,  then  the  two  shear  stress  curves 
shown  In  Fig.  1.6-2  would  coincide  exactly.  As  it  is,  the  degree  of 
agreement  attained  Is  considered  to  be  reasonably  satisfactory.  The 
discrepancy  between  the  two  curves  is  smaller  than  the  discrepancy  between 
the  experimental  points  and  either  of  the  curves. 

It  seems  probable  that  in  this  Instance  the  theoretically  computed 
values  are  actually  more  accurate  than  the  experimentally  measured  ones; 
the  experimental  measurement  Is  Innately  difficult  and  uncertain. 

It  may  be  concluded  that,  on  a  whole,. these  results  for  the  turbulent 

jet  substantiate  quite  well  the  proposed  heuristic  model  of  fluid  turbulence. 

* 

1.7  CONCLUSIONS  AND  RECOMMENDATIONS 

It  Is  concluded,  on  the  basis  of  the  evidence  available  so  far,  that 
the  proposed  heuristic  model  is  basically  adequate  for  determining  the 
principle  flow  characteristics  in  the  general  case  of  Inhomogeneous  and 
non-statlonary  turbulence  in  incompressible  flow. 

It  Is  recommended  that  the  present  model  be  applied  also  to  other 
cases  including  pipe  flow,  two  and  three  dimensional  wakes,  boundary  layers 
and  to  some  examples  of  unsteady  flow.  The  aim  would  be  to  refine  the 
present  unified  theory  and  extend  Its  range  of  applicability. 

It  Is  also  recommended  that  further  experimental  Information  be 
obtained  in  connection  with  those  aspects  of  the  model  for  which  the 
present  data  are  insufficient.  These  aspects  include,  for  example, 


(a)  the  relation  between  the  local  length  scale  X  of  the  mean  flow 
and  the  correlation  length  X*  of  the  turbulence, 

(b)  the  generalized  three  dimensional  stress/strain  rate  relations 
which  actually  exist  in  regions  of  strong  anisotropy  such  as 
in  the  flow  near  a  wall,  and 

(c)  the  influence  of  various  key  parameters  on  the  rate  of  dissipation 
of  turbulent  energy  into  heat. 
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T.  H.  Gawain,  D.Sc. 
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2.1  THE  FUNDAMENTAL  PROBLEM 


Consider  a  general  turbulent  flow  field,  one  which  Is  both  in¬ 
homogeneous  and  non-s tat ionary .  How  can  the  turbulent  structure  of 
such  a  field  be  adequately  analyzed  and  described? 

In  principle,  it  is  possible  to  formulate  a  valid  general  method 
for  computing  numerically  the  detailed  solution  u^n^(X,t)  for  a  particular 
realization  of  the  field,  subject  to  specified  boundary  conditions  and 
initial  conditions.  Once  these  conditions  are  fully  specified,  the 
equations  of  continuity  and  motion  suffice  to  establish  a  fully  determinate 
and  detailed  solution.  In  theory,  this  solution  procedure  could  be  repeated 
very  many  times  to  generate  a  whole  ensemble  of  such  detailed  realizations. 
All  realizations  would  be  compatible  with  the  same  general  macroscopic 
boundary  conditions  and  Initial  conditions.  However,  they  would  differ 
more  or  less  randomly  in  regard  to  the  fine  scale  details  of  the  boundary 
and  initial  conditions,  and  therefore  also  in  the  fine  scale  details  of  the 
resulting  motions. 

Once  such  an  ensemble  of  detailed  solutions  were  available,  any 
desired  statistical  features  of  the  turbulence  could  then  be  computed 
from  this  ensemble  in  a  completely  determinate  way. 

It  is  significant  that  the  foregoing  method  does,  in  principle, 
satisfy  the  mathematical  requirements  for  a  determinate  solution.  This 
solution  is  based  solely  and  6trictly  on  the  general  equations  of  contin¬ 
uity  and  motion  (the  Navier  Stokes  equations),  and  on  the  pertinent  bound¬ 
ary  and  Initial  conditions  for  the  particular  case  in  question. 
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In  practice,  however,  the  volume  of  detailed  calculations  required 
to  carry  out  the  above  hypothetical  scheme  of  calculation  for  any  case 
of  practical  Interest  is  simply  overwhelming.  It  lies  vastly  beyond  the 
capabilities  of  any  computer  yet  known  or  contemplated.  The  basic  reason 
for  this,  of  course,  is  that  the  detailed  space/time  structure  of  turbulence 
is  extremely  fine  eralned  and  enormously  intricate.  Hence  for  all  practical 
purposes,  the  above  hypothetical  solution  procedure  must  be  ruled  out. 

There  is  a  classical  method  for  attempting  to  circumvent  the  above 
difficulty.  Since  it  is  impractical  to  solve  the  detailed  equations  and 
average  the  results,  the  classical  method  averages  the  basic  equations 
themselves.  It  then  attempts  to  compute  the  desired  average  characteristics 
of  the  flow  from  these  averaged  equations.  Now,  if  the  basic  equations 
were  linear,  this  technique  would  indeed  solve  the  problem.  Unfortunately, 
they  are  not.  As  a  result  of  the  non-linearity,  the  averaging  process 
Introduces  additional  unknowns,  the  Reynolds  stresses,  into  the  basic 
equations.  Consequently,  there  are  finally  more  unknowns  than  available 
equations.  Therefore,  the  averaged  equations  do  not  comprise  a  closed 
and  determinate  set.  This  constitutes  the  so-called  closure  problem  of 
turbulence  theory. 

It  should  be  emphasized  that  the  deficiency  which  arises  from 
averaging  the  equations  is  fundamental.  No  amount  of  mere  manipulation 
of  the  equations,  however  Involved,  can  remedy  this  defect.  The  averaging 
of  the  equations  leads  to  a  definite  loss  of  essential  information  in 
comparison  with  the  information  implicit  in  the  original  equations. 

Therefore,  so  long  as  we  are  constrained  to  work  only  with  the 
averaged  equations  instead  of  with  the  originals,  it  becomes  unavoidably 
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necessary  to  Invoke  additional  hypotheses  in  order  to  compensate  for  the 
lost  Information  and  to  reestablish  a  determinate  set  of  equations. 
Unfortunately,  there  Is  no  known  way  to  establish  the  required  hypotheses 
deductively  from  the  original  fundamental  equations  themselves.  Consequently, 
It  Is  Initially  far  from  clear  just  what  auxiliary  hypotheses  are  most 
nearly  adequate  and  satisfactory. 

Obviously,  the  basic  justification  for  any  particular  hypotheses 
finally  adopted  must  therefore  rest  on  a  comparison  of  the  consequences 
deducts  from  these  hypotheses  with  experimental  observations.  This  is 
the  accepted  scientific  method  of  Inductive  inference.  While  it  Is  to  be 
regretted  that  the  problem  apparently  cannot  be  solved  on  a  wholly 
deductive  basis,  there  are  no  a-priorl  reasons  for  despairing  over  the 
prospects  for  ultimate  success  by  the  inductive  method. 

It  should  also  be  pointed  out  at  this  juncture  that  the  averaged 
equations  of  motion  describe  exactly  the  effect  of  the  Reynolds  stresses 
on  the  mean  motion.  However,  they  give  no  Information  whatever  concerning 
the  reciprocal  effect  of  the  mean  motion  on  the  Reynolds  stresses.  Yet 
some  definite  reciprocal  effect  of  this  kind  unquestionably  does  exist. 

This,  therefore,  is  precisely  the  problem  with  which  our  auxiliary 
hypotheses  must  necessarily  deal. 


2.2  STATISTICAL  DESCRIPTION  OF  TURBULENCE 


Consider  an  arbitrary  and  general  flow  field  which  is  neither 
homogeneous  in  space  nor  stationary  in  time  in  regard  to  its  turbulent 
structure.  Let  u£n\x,t)  denote  any  one  particular  realization  of  this 
field.  This  means  that  the  velocity  components  u^n^(Y, t)  represent  one 
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particular  detailed  solution  of  the  equations  of  continuity  and  motion; 
this  solution  Is  compatible  with  the  specified  macroscopic  boundary 
conditions  and  initial  conditions.  There  exist,  of  course  an  infinity  of 
such  possible  individual  realizations  which  satisfy  these  general 
conditions.  The  separate  realizations  differ  from  one  another  only  with 
respect  to  their  fine  scale  details.  We  denote  each  realization  by  a 
distinct  value  of  the  index  n,  which  is  taken  to  be  an  integer.  Thus  we 
may  have,  for  example,  n  ■  1,  2,  3,...N  where  N  is  some  indefinitely 
large  number.  Such  a  very  large  and  representative  collection  of  individual 
realizations  we  term  an  ensemble. 

Fortunately,  the  particular  details  of  any  one  individual  realization 
are  seldom  of  direct  interest.  What  is  usually  required  is  statistical 
information  concerning  the  principal  features  of  the  ensemble  as  a  whole. 

In  particular,  with  respect  to  any  variable  of  interest,  what  we  would 
generally  like  to  know  are  not  the  individual  values  for  particular 
realizations,  but  rather  the  following  items  of  information.  These  are 
listed  in  order  of  decreasing  urgency.  Thus 

(a)  The  average  value  of  the  variable  over  the  ensemble. 

(b)  The  statistical  variance  of  the  quantity  with  respect  to  the 
ensemble  average. 

(c)  The  detailed  probability  distribution  of  the  quantity  over 
the  ensemble. 

For  example,  if  Uj^n^Ot,t)  is  the  instantaneous  total  velocity  at 
space/time  point  x,t  for  realization  n,  we  require  the  corresponding 
ensemble  mean 

U.(x,t)-Umi£  U!(n)(x,t)  (2.2-1 

1  N-**>  n-1  1 
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This  result  defines  the  mean  flow  field.  Of  course,  the  turbulent 
fluctuations  of  the  nth  realization  are  defined  as  the  deviations  from 
this  mean .  Thus 

uJn)(K,t)  -  Uj(n)(JT,t)  -  U^X.t)  (-2) 

From  this  definition,  it  follows  that  the  turbulent  fluctuation  u^n^  has 
zero  mean  over  the  ensemble. 

The  corresponding  variance  of  this  velocity  fluctuation  is  given  by 

— 2 -  .  N  .  .2 

ui(x,t)  ■  iim  Tflljy  I  u^n  (x,t)  (No  sum  on  Index  i)  (-3) 

N-*»  n-1 

Finally,  the  probability  function  which  defines  the  detailed  statistical 
distribution  of  this  random  variable  u^  over  the  ensemble  may  be  denoted 
by  the  functional  4'[u^  (x,  t)  ] .  Here  the  quantity  Ydu^  represents  the 
probability  that  for  a  realization  n,  drawn  at  random  from  the  ensemble, 
the  velocity  component  u^  at  point  X,t  will  have  a  value  which  lies 
between  u^  and  (u^  +  du^) .  If  the  space/time  coordinates  x,t  be  hald 
constant  at  prescribed  values,  ¥(uj)  reduces  to  a  corresponding  ordinary 
probability  density  distribution.  Hence  the  functional  ¥  amounts  to  an 
ensemble  of  probability  distributions,  one  for  each  space/time  point  of 
the  field. 

If  the  turbulence  happens  to  b?  stationary  with  respect  to  time  in 
its  statistical  features,  then  it  becomes  permissible  and  convenient  to 
replace  the  ensemble  average  by  a  simple  time  average  over  a  single 
realization.  The  realization  index  n  may  therefore  be  dropped.  For 
example,  the  mean  velocity  distribution  may  now  be  redlflned  as  follows. 

1  +T 

U  (X)  -  Um -5=  /  u!(X,t)dt  (-4) 

1  T-k»  41  -1 
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Similar  changes  apply  to  the  definitions  of  variance  and  probability 
distribution. 

Likewise,  if  the  turbulence  is  homogeneous  with  respect  to  some 
line,  surface  or  volume,  the  ensemble  average  may  be  replaced  by  the 
corresponding  spatial  average. 

In  the  following  analysis,  the  average  will  always  be  an  ensemble 
average,  unless  specified  otherwise.  Also  the  discussion  will  center 
largely  on  certain  important  averages  and  variances,  hut  will  not  deal 
specifically  with  the  corresponding  detailed  probability  distributions. 


2.3  CORRELATION  FUNCTIONS 


In  order  to  formulate  an  adequate  theory  of  turbulence,  it  is 
necessary  to  define  certain  basic  concepts  which  characterize  the 
turbulence  in  the  general  vicinity  of  any  prescribed  space/time  point. 

In  this  connection,  the  two  point  velocity  correlation  tensor  is  of 
basic  importance. 

Let  P  denote  a  reference  space/time  point  with  coordinates  x,t 
regarded  as  arbitrary  but  fixed.  Let  P'  and  P"  be  two  associated  points 
with  coordinates 

x '  ■  x  +  Ax 
t'  -  t  +  At 

and  (2.3-1) 

x"  ■  x  -  Ax 
t"  -  t  -  At 

The  quantities  2Ax,  2At  represent  the  space  and  time  separation  between 
the  points  P'  and  P"  and  are  regarded  as  variable.  Note,  however  that 
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P*  and  P"  are  always  disposed  symmetrically  equal  and  opposite  with 
respect  to  the  fixed  reference  point  P. 

Components  of  the  turbulent  velocity  fluctuations  at  the  variable 
points  P'  and  P"  are  denoted  by 

uj^  »  u^(X  +  AX,  t  +  At)  at  P* 

and  (-2) 

uj  -  u^x  -  Ax,  t  -  At)  at  P" 

The  mean  product  of  these  two  velocity  components  is  the  tensor 

uJuT  -  uA(X  +  AX,  t  +  At)  Uj  (X  -  AX,  t  -  At)  (-3) 

where  the  overbar  denotes  an  ensemble  average. 

At  zero  separation,  the  above  tensor  reduces  (upon  reversal  of 
algebraic  sign)  to  the  ordinary  Peynolds  kinematic  stress  at  the  reference 
point  P.  Thus 

«  “^Uj  -  -uk(X,t)Uj(X,t)  (-4) 

By  reason  of  its  symmetry  the  stress  tensor  u^u^  =  u^u^  has  at  moat 
six  distinct  components.  The  general  velocity  correlation  tensor  u^u^', 
however,  is  not  necessarily  symmetrical  with  respect  to  an  interchange 
of  the  Indices  and  may  therefore  have  up  to  nine  distinct  components. 

Upon  contracting  the  Reynolds  stress  tensor,  we  obtain  the 
physically  important  quantity 

q2(x,t)  -  u1ui  -  u2  +  u2  +  u^  ■  2E(X,t)  (-5) 


Here  E  represents  the  ensemble  mean  turbulent  energy  at  the  point  x,t  and 
2 

q  is  t*e  corresponding  mean  square  velocity. 

We  now  define  the  space/time  correlation  tensor  R^j  merely  by 
normalizing  uju|j  with  respect  to  the  value  of  q2  at  the  fixed  reference 
point  P.  Thus 
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Pij  “  Rij^Jf»t5AX»At^ 


uju^  u1(X  +  AX,  t  +  At)u^ (X  -  AX,  t  -  At) 

2  2 
q  q  (5,t) 


It  Is  of  interest  to  compare  the  above  definition  which  we  hove 
adopted  with  certain  alternative  forms  which  are  sometimes  encountered. 
One  such  alternative  is 


(-6) 


R'  - 

1J 


u^(X,t)Uj(X  +  AX,  t  +  At) 


q(X, t)q(X  +  AX,  t  +  At) 
Another  alternative  is 


(-7) 


u  (X,t)u,(X  +  AX,  t  +  At) 

'/— — /—  - - 

y  ^(*,0 y Uj(x  +  ax,  t  +  At) 


(No  sum  on  1  or  j) 


(-8) 


The  form  we  have  chosen  to  preferable  because  it  is  related  in  a 
simpler  fashion  to  the  energy  spectrum  of  turbulence,  and  because  it  has 
certain  useful  symmetry  properties  which  will  be  pointed  out  presently. 

Two  Important  special  cases  of  the  foregoing  are  the. spatial  and 
temporal  correlation  tensors,  defined  respectively  as  follows: 


-  Ri:J(X,t;AX,0)  - 


R^  -  R^Oc.tjO.Ax) 


Uj(X  +  AX,  t)Uj(X  -  AX,  t) 


q2(X,0 


u1(X,t  +  At)Uj(X,t  ~  At) 
q2(S,t) 


(-9) 


Of  particular  significance  are  the  first  invariants  of  these  two 
tensors,  defined  as  follows: 


51 


(-10) 


,<•>  - 


R^“'  -  R^X.tjAx.O) 


ui(Jt  +  AX,  t)u^(X  -  AX,  t) 


q2(x,t) 


,(T)  _ 


Ril'  "  ^(*.^5°.^) 


Ui(X,t  +  At)u^(X,t  -  At) 


q2(*.t) 


Of  course,  the  summation  convention  applies  as  usual  to  the  Repeated 

indices  In  the  foregoing  expression.  Note  that  R^  Is  an  even  function  with 

(T) 

respect  to  the  Independent  variable  Ax.  Likewise  R^  is  even  with  respect  to 
At.  These  Important  symmetry  properties  are  a  consequence  of  defining 
Rjj  according  to  Eq.  (-6),  rather  than  according  to  the  alternative  (-7) 
or  (-8) . 

Because  of  their  fundamental  character  and  intimate  relation  to 

(b)  (T) 

the  turbulent  energy,  the  scalar  correlation  functions  R^  and  Rii 
provide  the  natural  bases  for  establishing  appropriate  length  and  time 
scales  of  the  turbulence  In  the  neighborhood  of  the  space/ time  point  X,t. 

Specifically,  let  dv'  ■  dxjdx^dx^  represent  an  infinitesimal 
element  of  volume,  and  dt*  an  Infinitesimal  element  of  time  at  the  variable 
point  P ' .  Also  let 

AS*  A*  -  AXjAXj^  ■  (AXj)2  +  (AXg) 2  +  (AXj) 2 
represent  the  square  of  the  spatial  separation  PP'. 

Then  the  desired  characteristic  length  X*  and  time  x*  may  be  defined 
as  follows: 


A* 


A*2(x,t) 


/|Rjj)|(AX*AX)dv' 


11 

/l^ldv' 


T*‘  - 


x*2(x,t) 


/|R 


™  I  (At)2dt' 


/IR±±> Idt' 


(-11) 


(-12) 
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The  integrals  In  (-11)  extend  over  all  space  at  a  fixed  time  t. 

The  integrals  in  (-12)  extend  over  all  tine  at  a  fixed  spatial  point  x. 

Of  course,  these  integrals  are  finite,  despite  the  infinite  domains  of 
integration,  because  of  the  rapid  decay  in  the  correlation  functions  at 
large  separations. 

Note  that  the  Integrands  contain  only  the  absolute  values  of  the 

correlation  functions;  this  ensures  that  every  non-zero  element  of  the 

correlation,  even  if  locally  negative  in  some  places,  makes  a  positive 

2  2 

contribution  to  the  Integrals  which  deff  t  X*  and  t*  . 

The  numerators  of  Eqs.  (-11)  and  (-12)  represent  the  second  moments 

of  the  respective  correlation  functions.  They  are  therefore  analogous 

2  2 

to  moments  of  inertia.  Also,  the  quantities  X*  and  t*  are  analogous  to 

variances  in  probability  theory.  It  is  well  known  that  such  second 

moments  attain  their  minimum  values  when  evaluated  with  respect  to  the 

centroldal  point  of  the  distribution  in  question.  In  the  present  instance 

the  correlation  functions  are  symmetrical,  with  their  centroids  at  the 

respective  origins  of  coordinates,  namely  at  Ax  »  0  and  at  At  ■  0.  Hence 

2  2 

the  characteristic  quantities  X*  and  t*  represent  the  minimum  possible 
"variances"  of  the  respective  correlation  functions.  These  properties 
make  X*  and  t*  particularly  appropriate  measures  of  the  length  and 
time  scales  of  turbulence  in  the  vicinity  of  the  arbitrary  space/time 
point  x,t. 


2.4  SPECTRUM  OF  TURBULENCE 

The  space/time  correlation  tensor  may  be  transformed  into  the 
wave  number/frequency  domain  by  means  of  multi-dimensional  Fourier 


53 


✓ 


integrals.  Let  tq  k 2  <3  be  the  components  of  the  spatial  wave  number  vec¬ 
tor  k.  Let  (o  be  temporal  wave  number,  that  is,  angular  frequency.  Also 
let  dv^  -  dtCjdtfjdtc^  represent  an  infinitesimal  element  of  "volume"  in 
wave  number  space.  Note  therefore  that 

R*AS  ■  -  tc^Ax^  +  +  <2^*3  (2.4-1) 

We  may  now  define  the  Fourier  transform  pair 

*i:j(*,t;*,u))  -  jj  Ri;j<x,tiAx,At)e+i(*‘A*  +  wAt)dv‘dt;  (-2) 

R^OqtjAS.At)  -  |j  •ljC*,tjlf,«)e~1(?,A*  +  “At)dvKdu  (-3) 

The  integrations  in  Eq.  (-2)  extend  over  all  of  space/time  A5c,At 
with  parameters  i,t  held  constant.  The  integrations  in  (-3)  extend  over 
all  of  wave  space  k,u  with  parameters  x,t  again  held  constant. 

The  spectral  correlation  function  4^  like  its  Fourier  transform 
R  is  a  nine  component  tensor.  Naturally  these  two  related  functions 
exhibit  corresponding  symmetry  properties. 

The  first  invarient  like  its  Fourier  transform  R^  is  a  scalar 
and  has  a  corresponding  basic  significance.  It  defines  the  character  of 
the  energy  spectrum  in  the  vicinity  of  the  space  time  point  x,t. 
Incidentally,  it  can  be  shown  from  energy  considerations  that  must  be 
non-nega  t ive  everywhere . 

It  may  be  seen  that  the  functions  and  9^  are  exactly  equivalent 
in  the  information  they  provide  concerning  the  structure  of  the  turbulent 
field.  The  reason  is  simply  that  if  either  one  of  these  functions  be 
specified,  the  other  one  can  be  found  from  the  above  Fourier  Integral 
relations. 
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It  Is  clear  therefore  that  fundamental  scales  of  length  and  time, 
which  were  previously  developed  in  terms  of  R^®'  and  R^;  can  equally 
well  be  defined  in  terms  of  the  corresponding  transforms  and  . 

Because  of  the  above  equivalence  between  the  correlation  functions 
Rjj  and  the  spectral  functions  it  is  immaterial  in  terms  of  which  of 
these  our  turbulence  theory  is  developed.  We  have  chosen  to  emphasize 
the  correlation  aspects,  and  it  is  seen  that  there  is  no  loss  of  infor¬ 
mation  Involved  in  so  doing. 

It  should  be  noted  that  the  nine  spectral  functions  F^OCftjlX.At) 
are  ostensibly  eight  dimensional,  the  Independent  variables  being  x^,  *  , 
x^,  t,  Ax^,  AXj,  Ax^,  At.  Similarly,  the  nine  spectral  functions  4^ 
(R,t;R,w)  are  also  eight  dimensional,  the  independent  variables  being 

Xj.  *2*  x3»  *^1*  k2"  *3’ 

In  addition  to  the  above  eight  Independent  dimensions  which  are 
explicit,  there  is  an  additional  dimension  which  is  implicit.  Recall 
that  the  correlation  functions  are  themselves  averages  theoretically 
computable  from  a  large  ensemble  of  individual  flow  realizations.  We  may 
assign  to  each  realization  a  superscript  index  n  which  varies  from  one  to 
some  very  large  number.  Then  n  is  equivalent  to  an  additional  Independent 
dimension.  For  if  we  wish  to  know  not  just  the  ensemble  averages  them¬ 
selves,  but  also  something  about  their  statistical  distributions,  then  it 
becomes  necessary  to  take  into  account  the  possible  variations  of  the 
Independent  realization  n  from  the  ensemble  mean. 

On  the  other  hand,  any  specialization  of  the  conditions  reduces 
the  number  of  Independent  variables  which  must  be  considered.  Thus,  for 
stationary  turbulence,  R^  becomes  Independent  of  t  (but  not  of  At)  For 
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completely  homogeneous  turbulence,  R^  becomes  Independent  of  x^,  x^,  x^. 

Also,  for  Isotropic  conditions,  spherical  symmetry  prevails,  and  R^ 
ceases  to  depend  on  Ax^,  AXj,  Ax^  Individually,  but  becomes  a  function 
only  of  their  resultant  ^Ax^) 2  +  (AXj)^  +  (Ax^)^ •  Furthermore,  in  this 
case  the  number  of  distinct  components  of  reduces  from  nine  to  two, 
depending  only  on  whether  the  Indices  l,j  are  eoual  or  unequal.  These 
are  vast  simplifications.  They  account  for  the  fact  that  homogeneous 
and  Isotropic  turbulence  has  been  studied  so  much  more  extensively 
than  any  other  case. 

It  should  be  noted  that  the  fundamental  phenomena  of  turbulence  aa 
represented  by  the  ensemble  of  functions  u^(x,t)have  five  Independent 
variables,  namely,  x^,  x«,  t  and  n,  and  three  dependent  variables,  namely, 

u^,  Uj,  and  u^.  It  Is  a  curious  fact  that  when  we  seek  to  make  these 
phenomena  more  comprehensible  by  means  of  the  correlation  functions  R4 
the  number  of  dependent  and  Independent  variables,  instead  of  being 
sharply  reduced,  are  actually  Increased!  Of  course,  the  functions  R^ 
have  a  regular  and  symmetrical  structure,  they  vary  smoothly,  and  die 
away  rapidly  at  large  separations,  whereas  the  u^  vary  in  an  extremely 
complex  and  erratic  way.  It  is  this  qualitative  difference  that  makes  the 
correlation  functions  seem  so  much  more  comprehensible. 


ij 


2.5  RELATION  BETWEEN  TURBULENT  CORRELATION  SCALES  AND  MEAN  FLOW 


It  Is  customary  and  in  many  ways  convenient  to  resolve  fluid 
turbulence  Into  two  distinct  complementary  aspects  which  we  then  refer  to 
as  "the  mean  flow"  and  "the  turbulent  fluctuations",  respectively.  While 
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such  a  distinction  Is  useful  for  certain  conceptual  and  computational 
purposes.  It  should  not  be  forgotten  that  these  two  aspects  are  inextri¬ 
cably  coupled  within  a  single  unified  process.  In  this  connection  it  Is 
pertinent  to  note  that  while  we  often  find  it  convenient  to  introduce 
this  distinction  between  mean  flow  and  turbulent  fluctuation  into  the  basic 
Navler  Stokes  equations,  the  distinction  is  not  Innate  within  the  original 
equations  themselves.  This  Is  shown  by  the  fact  that  these  equations, 
being  non-linear,  cannot  be  fully  decoupled  into  two  Independent  sets, 
one  dealing  exclusively  with  the  mean  flow,  the  other  dealing  exclusively 
with  the  turbulent  fluctuations. 

The  foregoing  fact  provides  a  useful  hint  In  connection  with  the 

problem  of  formulating  appropriate  auxiliary  hypotheses,  as  required  for 

supplementing  the  averaged  equations  of  motion.  We  have  seen  that  any 

such  hypotheses  will  naturally  be  formulated  in  terms  of  the  fundamental 

physical  parameters  which  characterise  the  overall  turbulence  In  the 

vicinity  of  a  given  space/time  point  5c , t .  Among  these  quantities  are 

2 

the  mean  square  velocity  q  (x,t)  and  the  macroscales  of  turbulence  in 
space  and  time,  as  expressed  by  the  previously  defined  parameters  A*  and 

T*. 

Of  course,  In  most  problems  of  practical  Importance,  Interest  centers 
chiefly  on  the  mean  flow  field.  One  reason  for  this  is  that  the  mean 
flow  Is  usually  much  more  evident  physically,  much  more  convenient  to 
observe  and  measure,  seemingly  more  easy  to  understand.  We  tend  to 
regard  the  turbulent  fluctuations  as  a  nuisur.ce  which  has  to  be  dealt  with 
In  order  to  get  to  the  mean  flow,  and  we  tend  to  confine  attention  to 
those  aspects  of  the  turbulence  whose  connection  to  the  mean  flow  is 

V 

obvious  and  unavoidable. 
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More  often  than  not,  we  do  not  have  available  the  data  which  would 
be  needed  to  establish  such  largely  unseen  characteristics  of  the  turbul¬ 
ence  as  the  spatial  macroscale  X*  and  temporal  macroscale  t*.  Therefore, 
it  is  of  doubtful  utility  to  formulate  any  heuristic  theory  of  turbulence 
directly  in  terms  of  these  quantities.  For  example,  what  good  does  it 
do  to  express  the  Reynolds  stresses,  say,  in  terms  of  X*  and  t*  if  the 
latter  quantities  are  initially  just  as  much  unknown  as  the  Reynold's 
stresses  themselves?  Consequently,  we  must  at  this  stage  of  our  knowledge 
replace  the  fundamental  macroscales  X*  and  t*  in  our  formulation  by  other 
quantities  which  are  more  or  less  equivalent,  but  which  are  far  more 
accessible  to  observation  and  calculation. 

Fortunately,  for  this  purpose  we  can  take  advantage  of  the  intimate 
connection  as  noted  above,  which  necessarily  exists  between  the  local 
turbulence  field  and  the  local  mean  flow  field.  Presumably  because  of 
this  connection,  these  local  scales  of  space  and  time  must  somehow  be 
reflected  in  the  local  characteristics  of  the  mean  flow,  and  must 
therefore  be  deduclble  from  the  mean  flow. 

While  the  exact  nature  of  this  relationship  between  the  macroscales 
and  the  mean  flow  is  far  from  clear  we  can  legitimately  infer  that  the 
relationship  is  inherently  somewhat  local  in  character.  This  follows  from 
a  fact  universally  revealed  by  all  observed  correlation  curves.  These 
always  show  a  rapid  decrease  in  correlation  with  increasing  separation  in 
space  and  time.  Moreover,  the  correlations  become  negligible  at  separa¬ 
tions  larger  than  about  two  or  three  multiples  of  X*  or  r*. 

Consequently,  we  infer,  or  at  least  we  hypothesize,  that  quantities 
more  or  lesc  equivalent  to  the  true  macro  scales  X*  and  t*  can  be  defined 
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and  computed  In  terms  of  the  observable  characteristics  of  the  local  mean 
flow.  Let  us  call  these  quantities  the  apparent  macroscales  and  designate 
them  by  the  symbols  A  and  t  to  distinguish  them  from  the  true  macroscales 
A*  and  t*  as  earlier  defined.  The  Important  problem  of  just  hov  A  and  t 
should  be  defined  in  terms  of  the  observable  mean  flow  characteristics 
is  taken  up  In  a  subsequent  section.  For  the  present  we  observe  simply 
that  if  A  and  t  be  defined  in  some  suitable  manner,  then  it  seems  probable 
from  the  above  reasoning  that  the  dimensionless  ratios 


—  B  f 

A*  A 


(2.5-1) 


7*  -  f T  <“2> 

will  be  very  stable  and  well  behaved  universal  functions  which  vary  but 
slowly  if  at  all.  We  assume  that  for  many  practical  applications  it  will 
turn  out  to  be  premlsslble  to  approximate  these  universal  functions  by 
suitable  average  values,  and  to  treat  them  as  simple  quasi-constants.  The 
actual  numerical  values  of  these  ratios  is  regarded  as  a  matter  to  be 
determined  empirically.  Of  course,  these  numerical  values  will  depend  on 
the  exact  definitions  that  are  finally  adopted  for  A  and  t.  However,  the 
presumption  is  that  if  the  definitions  adopted  are  near  optimal,  the 
ratios  will  turn  out  to  be  roughly  of  order  unity. 

It  follows  from  the  foregoing  discussion  that  our  definitions  of  the 
apparent  macroscales  A(x,t),  i(x,t)  in  the  vicinity  of  an  arbitrary  space/ 
time  point  x,t  should  depend  only  on  mean  flow  conditions  in  a  finite 
region  surrounding  that  point.  This  principle  is  more  realistic  than  any 
method  which  seeks  to  define  these  macroscales  in  terms  of  mean  flow 
quantities  only  at  the  point  x,t  itself.  On  the  other  hand,  it  avoids 
the  opposite  extreme  which  would  hold  that  A(S,t)  and  T(S,t)  are  somehow 
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directly  dependent  on  conditions  at  all  points  in  the  flow  field,  however 
remote  in  space  or  time.  Our  proposed  model  accords  more  nearly  than 
either  of  these  alternatives  with  the  observed  fact  that  turbulent 
correlations  are  effectively  limited  to  certain  finite  separations. 

Naturally,  the  definitions  of  X(£,t)  and  i(x,t)  should  emphasize 
mean  flow  conditions  in  the  Immediate  vicinity  of  the  point  X,t,  should 
give  progressively  less  weight  to  conditions  farther  away,  and  should 
finally  neglect  mean  flow  effects  at  points  sufficiently  remote  in  space 
or  time  from  point  x,t.  This  requirement  suggests  the  general  idea  of 
defining  X  and  t  in  terms  of  appropriate  weighted  averages  of  certain 
pertinent  mean  flow  quantities. 

The  preferred  choice  of  weighting  functions  for  this  purpose  would 
be  the  correlation  functions  and  themselves,  if  these  were 

known.  Of  course.  If  these  functions  were  known,  they  could  be  used  to 
find  the  true  macroscales  X*  and  t*  directly,  and  there  would  then  be  no 
need  for  weighting  functions.  However,  they  are  not  known.  Consequently, 
we  are  obliged  to  used  heuristic  substitutes,  preferably  functions  which 
resemble  the  correlation  functions  to  a  certain  extent.  Fortunately, 
great  accuracy  is  not  necessary  in  this  regard,  for  weighted  averages 
tend  to  be  relatively  insensitive  to  minor  variations  in  the  form  of  the 
weighting  functions  used. 

In  view  of  these  considerations,  wc  choose  the  Gaussian  curve  as 
an  appropriate  form  of  heuristic  weighting  function.  It  expresses  the 
general  trends  of  correlation  in  a  suitable  manner,  and  has  various 
convenient  and  known  mathematical  properties  as  well.  Consequently,  we 
finally  write  the  space/time  weighting  function  as  follows 
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Ax»Ag  (At)  2 
2  2 

W  »  e  T  (-3) 

The  reason  for  inserting  the  numerical  factor  3  in  the  first  term  of  the 
exponent  will  be  explained  presently. 

Of  course,  for  purely  spatial  weighting,  we  may  delete  the  temporal 
term  in  the  exponent  and  set 


_3 

w<*>  -  .  »2 

Similarly,  for  purely  temporal  weighting,  v;e  may  delete  the  spatial  term 
in  the  exponent  and  set 


<-*) 


.  (^t)2 

W<*>  -  e 

(s)  (T) 

The  above  weighting  functions  Wv  1  and  VT  '  may  be  regarded  as 

(  g\ 

initial  heuristic  approximations  to  the  true  correlation  functions 
(T) 

and  ,  respectively.  Note  that  all  four  of  these  functions  equal 
unity  at  zero  separation. 

If  we  now  replace  the  true  correlation  functions  by  their  heuristic 
Gaussian  approximations  in  Eqs.  (2.3-11)  and  (2.3-12),  we  obtain  the 
interesting  results 


(-5) 


2  2 

A*  -  A 

(-6) 

*2  2 

T*  ■  T 

(-7) 

This  shows  that  for  Gaussian  correlation  functions  with  exponential 
constants  A  and  t,  the  true  macroscales  A*  and  t*  may  be  simply  identified 
with  these  exponential  constants.  Of  cwurse,  this  simple  identification 
is  possible  only  if  the  numerical  factor  3  is  inserted  into  the  spatial 
exponent  as  indicated.  This  arises  from  the  fact  that  the  spatial  integra¬ 
tion  of  Eq.  (2.3-11)  is  three  dimensional,  whereas  tha  temporal  integration 
of  Eq.  (2.3-12)  is  only  one  dimensional. 
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Recall,  however,  that  the  actual  correlation  functions  are  not 
necessarily  Gaussian  and  are  in  fact  unknowns  in  the  present  analysis. 
Moreover,  for  the  reasons  previously  given,  the  quantities  X  and  t  will 
be  defined  uot  in  terms  of  the  correlation  functions,  but  heurlstically 
in  terms  of  the  mean  flow.  In  fact  the  basic  reason  for  introducing  the 
above  weighting  functions  in  the  first  place  is  that  they  are  needed  in 
connection  with  these  heuristic  definitions  of  X  and  t.  However,  the 
weighting  functions  in  turn  contain  X  and  t  as  parameters.  It  is  there¬ 
fore  apparent  that  the  actual  calculation  of  X  and  t  will  necessarily 
involve  some  type  of  iteration  process. 

For  certain  purposes  in  the  subsequent  analysis  it  will  prove 
convenient  to  utilize  a  normalized  form  of  the  above  Gaussian  weighting 
function.  In  this  connection,  consider  the  following  integrals,  namely. 


H  -  //Wdv'dt'  (-8) 

H(s)  -  //W^dv'dt'  (-9) 


and 

H<T>  -  //W(T)dv'dt  (-10) 

As  usual,  the  integrands  extend  over  the  entire  space/ time  domain.  Of 
course,  the  integrals  themselves  are  always  finite  even  if  the  domain  of 
Integration  extends  to  infinity. 

We  now  define  the  normalized  weighting  functions  as  follows: 


w 


W 

H 


v 


(s) 


W<8) 

73 


w 


(T) 


(T) 


I! 


(-11) 

(-12) 


(-13) 
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When  normalized  in  this  way,  the  above  weighting  functions  have  the 
following  convenient  properties,  that  is, 

//wdv'dt'  ■  1  (-14) 

//w^dv'dt'  ■  1  (-15) 

//w(T)dv'dt’  -  1  (-16) 

2.6  LOCAL  LENGTH  AND  TIME  SCALES  OF  THE  MEAN  FLOW 

We  have  been  considering  any  arbitrary  and  general  turbulent  flow 
field  of  which  the  mean  flow  may  be  either  steady  or  unsteady.  With  every 
space/time  point  of  this  field  may  be  associated  two  parameters  X  and  t 
which  express  respectively  a  length  scale  and  a  time  scale  characteristic 
of  the  mean  flow  pattern  ir.  the  vicinity  of  the  point.  Our  problem  is  to 
devise  suitable  explicit  definitions  of  these  parameters  in  a  way  which 
meets  various  essential  requirements.  Some  of  these  have  already  been 
discussed.  There  are  some  additional  requirements  as  well,  among  which 
are  the  following: 

(1)  The  parameters  X  and  t  must  be  true  scalars,  and  therefore  lnvarlent 
with  respect  to  any  rotation,  reflection,  translation  or  acceleration 
of  the  reference  axes. 

(2)  These  parameters  should  also  be  everywhere  continuous,  finite  and 
positive  (except  possibly  in  certain  limited  singular  regions,  such 
as  at  a  solid  boundary) . 

2  2 

(3)  If  feasible,  the  mean  flow  quantities  X  and  t  .  should  preferably  be 

2  2 

related  to  the  correlation  quantities  X*  and  t*  ,  respectively,  such 
X2  2 

that  the  ratios  — r  and  — are  close  to  unity,  or  at  least  so  that 

X*^  T* 

these  ratios  remain  as  nearly  constant  as  possible. 
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These  requirements  can  be  satisfied  as  follows.  Let 

Ui  "  Ui^xi'x2,x3»t^  "  (2.6-1) 

represent  the  velocity  field  of  the  mean  flow.  Then  the  corresponding 
strain  rate  tensor  is 


r 


ij 


(-2) 


Of  course,  we  are  considering  here  only  the  incompressible  case, 
for  which  the  first  invariant  vanishes,  that  is 


(-3) 


Consequently,  there  is  no  net  dilatation,  and  the  1*^  represent  purely 
distortional  effects. 

We  now  define  a  generalised  strain  rate  0,  and  a  generalized  strain 
rate  gradient  O'  by  means  of  the  following  expressions,  namely 


lr  r 

TVij 


<-'♦) 


n 


,2 


(-5) 


From  these  definitions,  the  followinp  relation  can  be  obtained.  It 
will  prove  useful  in  the  subsequent  work.  Thus 


(-6) 


The  factor  1/2  which  appears  in  Eq.  (-4)  calls  for  a  brief  explanation. 
Consider  a  simple  plane  parallel  shear  flow  such  as  illustrated  in  Fig. 

2.6-1.  Let  the  velocity  be  parallel  to  axis  x^,  and  let  the  shear  gradient 
be  in  the  direction  of  x^.  Then  ■  0.  Now 

the  only  non-zero  shear  strain  rate  is  ■  r^.  Then  applying  Eq.  (-4) 


64 


y 


In  this  com  the  definitions  of  SI  and 
SI'  reduce  to 

*U  ■  (dU/dy) 

SI'*  U"  *  (d*U/dy*) 


FIGURE  2.6-1  EXAMPLE:  SHEAR  3UANTITIES  SI  AND  SI' 
IN  PLANE  RARALLEL  FLOW 
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we  find  that,  because  of  the  presence  of  the  factor  1/2,  we  obtain  the 
simple  result 


ft  -  r 


12 


In  other  words  the  presence  of  the  optional  factor  1/2  identifies  ft 
with  the  ordinary  shearing  rate  in  a  simple  plane  shear  flow.  Eq ,  (-4) 
may  be  regarded  as  an  extension  of  this  comcept  to  a  general  state  of 
strain  rate. 

It  can  be  shown  that  while  and  depend  on  the  orientation  of 

the  reference  axes,  the  quantity  Q  is  invarlent  in  this  regard,  and  is 

12  2 

therefore  a  true  scalar.  The  same  is  true  of  the  quantities  ft  ,0  ,  ft ' 

2  2  2 

and  (nfl1)  .  Hence  X  and  t  can  be  conveniently  defined  in  terms  of  any 
Independent  pair  of  these  quantities. 

To  illustrate  this  idea,  consider  the  expressions 


X 


2 


(-8) 


2  2  1 
T  -c’? 


2  2 

where  and  C  are  dimensionless  constants.  The  combinations  shown  may 
be  seen  to  define  local  scales  of  length  and  time,  X  and  t.  In  fact, 

Eq.  (-8)  will  be  recognized  as  the  expression  for  the  von  Karman  mixing 
length,  or  rather,  as  a  generalization  of  the  concept,  for  von  Karman's 
original  definition  was  restricted  to  the  special  case  of  parallel  flow. 

One  difficulty  with  the  above  approach,  however,  is  that  the  functions 

■j  2 

X  and  t  so  defined  tend  to  fluctuate  too  wildly  Instead  of  being  smooth 

and  well  behaved.  Moreover,  at  certain  singular  points  where  $  or 

2  2 

are  locally  zero,  the  quantities  X  or  t  may  reach  values  of  zero  or 
irfinity,  neither  of  which  is  prenissible  or  physically  meaningful. 
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The  cure  for  this  difficulty,  of  course,  is  to  replace  the  local 
2  2 

point  values  of  R  and  (1 1  by  suitable  space/time  averages  evaluated  in 
the  vicinity  of  the  point.  Thus  we  may  write 


R2  -  //wfi2dv'dt ' 
R*2  ■  //wR  ,2dv'dt ' 

whereupon 


(-10) 

(-11) 

(-12) 


2  2 

The  new  values  of  X  and  t  obtained  in  this  way  will  now  vary 

smoothly  as  required,  and  will  never  be  either  zero  cr  infinite.  The 

2  2 

dimensionless  constants  C.  and  C  remain  to  be  specified. 

A  T 

Note  that  the  normalized  form  of  the  weighting  function  is  used  in 
the  Integrals  (-10)  and  (-11) .  This  is  Important  because  it  ensures  t  lat 
the  Integrals  remain  relatively  insensitive  to  the  precise  form  of  the 
weighting  function. 

For  the  reasons  given  earlier,  a  Gaussian  curve  has  been  chosen  for 
the  weighting  function.  For  consistency,  the  characteristic  exponential 
constants  X  and  t  of  this  function  should  be  respectively  equal  to  the 
X  and  t  computed  from  F.qs.  (-12)  and  (-13).  In  other  words,  the 
averaging  intervals  in  length  and  time  should  be  approximately  equal  to 
the  characteristic  local  scales  of  length  and  time.  Some  iteration  will 
be  Involved  in  satisfying  this  requirement. 

It  turns  out  to  be  useful  to  generalize  the  above  idea  even  further. 
For  this  purpose  we  define  the  following  integrals,  namely, 


i2  -  //v(n2)a(fl,2)bdv'dt' 


(-14) 


J2  -  //«(fi2)C(^J,2)ddv,dt,  (-15) 

where  the  exponents  a,  b,  c,  d  may  be  assigned  arbitrarily.  However, 
for  reasons  which  will  soon  be  evident,  we  stipulate  that 


D 


b  d 

(a  +  b)  (c  *f  d) 


f  0 


(“IS) 


Moreover,  to  avoid  improper  integrals  we  require  all  of  the  above  exponents 
to  be  non-negative. 

How,  depending  only  on  how  the  arbitrary  exponents  a,  b,  c,  d  are 

2  2 

specified,  the  Integrals  1  and  J  may  be  assigned  various  physical 

dimensions.  Consequently,  it  is  always  possible  to  define  a  mean  square 
2  2 

length  X  and  mean  square  time  i  by  expressions  of  the  following  general 
form,  namely, 

A2  -  C2(I2)k(jV  (-17) 

t2  -  C2(I2)®(J2)"  (-18) 

2  2 

where  and  are  arbitrary  dimensionless  constants.  It  Ip  easily  shown 
on  dimensional  grounds  that  the  four  exponents  k,  A,  m,  n  are  given  by 
the  following  formulas: 


k  -  -  ~(c  +  d) 


l  -  +  ^(a  +  b) 


(-19) 


It  is  now  seen  that  Eq.  (-16)  represents  a  necessary  condition  to 
ensure  the  existence  of  the  above  solution. 
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The  foregoing  formulas  show  that,  subject  only  to  the  condition  (-16), 

2  2 

specification  of  the  six  arbitrary  constants  a,  b,  c,  d,  C.  and  C  uniquely 

A  T 

2  2 

defines  the  quantities  1  and  t  .  These  six  auantlties  are  here  regatded 

as  empirical  constants  whose  values  can  in  principle  be  chosen  so  as  to 

2  2 

confer  upon  the  defined  quantities  X  and  t  a  wide  range  of  desired 
properties. 

Ideally,  it  would  be  desirable  to  choose  these  six  constants  in 

2  2  2 
such  a  way  as  to  make  X  ant  t  agree  as  nearly  as  possible  with  X*  and 

2 

t*  ,  respectively,  over  a  broad  range  of  conditions.  In  principle,  this 
can  always  be  done  if  and  when  sufficient  experimental  data  becomes 
available  for  this  purpose.  At  the  present  time,  the  available  data  are 
Insufficient  for  such  an  analysis.  It  therefore  becomes  necessary,  at 
least  for  the  time  being,  to  resort  to  expedients  of  a  more  heuristic 
nature,  as  will  be  explained  presently. 

2  2  2 

The  calculation  of  the  space/time  Integrals  I  and  J  which  fix  X 

2 

and  t  is  seen  to  be  inherently  lengthy  and  complicated.  One  reason  is 

that  the  space/ time  integrals  involved  are,  in  their  present  form,  four 

dimensional.  The  second  reason  is  that  these  samt  quantities  appear  not 

only  as  the  final  unknowns,  but  also  as  the  exponential  constants  in  the 

Gaussian  weighting  function.  Hence  our  next  endeavor  is  to  simplify  the 

analysis  as  much  as  possible  in  both  these  respects. 

In  this  connection,  we  observe  first  that  in  computing  the  integrals 
2  2 

I  and  J  ,  if  the  mean  flow  happens  to  be  steady,  the  Integration  over 
time  cancels  identically  from  the  result,  and  the  integrals  reduce  to 
ordinary  volume  integrals.  Secondly,  we  notice  that  if  the  flow  field  also 
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happens  Co  be  either  plane  or  axl-symmetrlc,  all  volume  Integrals  can  be 
reduced  tc  corresponding  surface  Integrals.  Also  If  the  flow  field  happens 
to  be  of  the  type  known  as  self-similar,  all  volume  integrals  again 
reduce  to  surface  Integrals.  If  the  flow  field  is  both  ax  l-syrametric  and 
self-similar,  then  a  double  reduction  is  possible,  and  the  volume  integrals 
can  be  reduced  to  line  integrals.  These  are  tremendous  simplifications. 

By  a  self -similar  flow  field,  we  mean  a  field  in  which  there  exists 
a  family  of  surfaces  across  which  the  mean  velocity  distribution  has  a 
common  pattern;  the  velocity  distribution  across  the  various  surfaces 
differs,  at  most,  only  in  the  respective  scales  of  length  and  time,  or 
what  amounts  to  the  same  thing,  of  length  and  velocity.  Many  common  flows 
in  jets,  wakes,  duett-.:  and  boundary  layers  are  known  to  be  self-similar. 

The  characteristic  surfaces  in  such  cases  are  conveniently  taken  as  cross- 
sections  approximately  normal  to  the  principal  flow  direction. 

This  consideration  of  self-similarity  shows  us,  thirdly,  that  for 
such  flows  each  characteristic  surface  possesses  its  own  overall  length 
and  time  scales;  these  are  constant  over,  any  one  such  surface,  but  may 
differ  from  one  surface  to  the  next.  In  addition  to  these  constant 
overall  scales  of  length  and  time  which  characterize  the  cross-section  as 
a  whole,  there  arc  also  distinct  local  scales  of  length  and  time,  as 
considered  earlier,  in  the  vicinity  of  any  reference  point  on  the  surface. 
These  local  scales,  unlike  the  overall  scales,  vary  smoothly  from  point  to 
point  along  the  cross-sectional  surface. 

Of  course,  even  when  we  have  taken  these  factors  into  account  to 

2  2 

simplify  our  formulas,  the  definitions  of  X  and  t  will  remain  indeter¬ 
minate  until  the  four  exponents  a,  b,  c,  d  and  the  two  proportionality 
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2  2 

constant*  C£  and  are  definitely  specified.  Obviously,  there  exists  an 
enormous  amount  of  flexibility  here,  and  this  was  done  Intentionally. 

The  flexibility  has  been  deliberately  provided  in  order  to  permit  matching 

2  2  2  2 
X  and  t  as  closely  as  possible  to  X*  and  t*  over  a  broad  range  of 

conditions. 

Naturally,  to  do  this  properly  will  require  much  data  beyond  what  is 
now  available.  Hence  at  the  present  time  our  scheme  has  far  more 
flexibility  than  we  really  need  or  can  properly  exploit. 

Consequently,  we  have  found  it  necessary  to  simplify  at  this  time  in 
what  may  seem  to  be  a  somewhat  arbitrary  manner.  Actually,  such  simplifi¬ 
cations  were  guided  by  a  fair  amount  of  numerical  experimentation.  Accord¬ 
ing  to  present  evidence,  they  seem  to  give  reasonably  satisfactory  agree¬ 
ment  with  physical  experiment.  No  doubt,  there  remains  room  for  further 
improvement,  however. 

Specifically,  for  our  present  turbulence  model,  we  have  chosen  the 
following  numerical  values  for  the  four  arbitrary  exponents,  namely, 

a  ■  2  b-0  c-1  d  -  1  (-20) 

Also,  it  proves  to  be  most  convenient  to  choose  the  proportionality 
‘ factor  simply  as 

Cx2  -  1  (-21) 

2 

The  other  proportionality  factor  is  not  required  for  the  present  simpli¬ 
fied  analysis. 

Let  us  summarize  the  situation  as  it  now  stands  after  the  foregoing 
simplifications.  The  integrals  below  extend  over  the  entire  flow  field. 

Thus 

Lit  •  (it'  -  <)  -  separation  variable  (-25) 
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W(*,*')  -  e 


weighting  function 


Lit' Lit 


H(<)  ■  / Wdv'  ■  normalizing  integral 
v(t,it)  ■  ^  W  ■  normalized  weighting  function 

I2(<)  ■  /wfl^dv'  ■  first  characteristic  integral 
2  2 

J  (it)  ■  /w(nn')  dv'  ■  second  characteristic  integral 
2  I2 

A  (<)  ■  — «  ■  local  length  scale  parameter 


(-26) 

(-27) 

(-28) 

(-29) 

(-30) 

(-31) 


Notice  that  Eq.  (-26)  above  differs  from  Eo.  (2.5-4)  in  that  the 
unessential  factor  3  has  been  eliminated  from  the  exponent  for  convenience. 
Thitt  is  arbitrary  but  permissible;  it  amounts  to  a  rescaling  of  A. 

The  foregoing  relations  suffice  in  principle  to  define  the  required 
scale  factor  A  at  every  point  in  the  flow  field  as  required,  but  do  involve 
a  computational  difficulty.  The  problem  is  that  the  value  of  A  is  needed 
for  Eq.  (-26)  before  a  computed  value  is  available  from  Eo.  (-31).  Hence, 
it  is  necessary  to  resort  to  an  iteration  process  such  as  the  following. 

Let  the  (n  +  l)th  approximation  to  A  be  defined  as: 


»Ll®  ■  *Ll®"Ll® 

(-32) 

where 

'Ll®  ■  Un 

. (-33) 

JLi®  •  /. n 

(-34) 

The  weighting  function  is  defined  as  follows.  Let 

W  (5f,5*f)  -  exp  -  (*'  ~ 

(-35) 
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Then 


_  +  W  (x,x) 

w  (x,x)  ■  - 

f  ..  W  (x,x')dv' 

Jall  space  n  * 


(-36) 


We  now  define  (in  principle) 


A^(x)  ■  tim  A^(x) 
n~  n 


(-37) 


Advantage  should  be  taken,  of  course,  of  any  knowledge  concerning  the 
flow  to  hasten  the  convergence  of  the  A^'s.  In  particular,  if  the  flow  is 
self-similar,  such  as  the  flow  in  a  pipe,  duct,  or  boundary  layer,  or  the 
far  downstream  region  of  a  turbulent  jet  or  wake,  calculations  need  only 
be  performed  on  a  single,  representative  cross-section  of  the  flow,  and 
the  results  rescaled  according  to  the  appropriate  nondlmenslonallzlng 
parameters. 

As  an  initial  guess  to  start  the  iteration  process,  we  take  Aq  ■  “ 
everywhere.  Then  in  computing  A ^ ,  we  find  that  the  weighting  function 
Wq  is  equal  to  unity  everywhere,  and  so  we  obtain  simply 


x2  „  all.  spacg -  (-38) 

/ail  space  dv 

a  constant  independent  of  position.  Thus,  is  the  first  nonconstant 

approximation  to  A  that  we  obtain. 

Numerical  experiments  with  various  flow  fields  have  shown  that,  in 

general,  the  convergence  of  the  A  's  is  extremely  rapid.  This  is  illustrated 

n 

in  Table  2.6-1  which  shows  successive  numerical  approximations  to  the  A 
distribution  in  laminar  flow  in  a  two  dimensional  channel.  Consequently, 
for  the  purposes  of  the  calculations  presented  in  this  paper,  the 
approximation  was  made: 


A(x)  s*  A2(x) 


(-39) 


Some  improvement  could  possibly  be  gained  by  continuing  to  higher 
approximations,  but  in  view  of  the  gross  simplifications  already  inherent 
in  the  model  as  a  whole,  the  authors  feel  that  such  refinement  is  not 
warranted.  Moreover,  absolute  accuracy  in  this  regard  is  not  essential. 
It  is  necessary,  however,  that  the  chosen  order  of  approximation  be 
adhered  to  consistently.  It  is  also  important  to  adjust  the  initially 
undetermined  numerical  constants  of  the  model  so  that  the  chosen  order  of 
approximation  gives  satisfactory  agreement  with  experimental  data.  The 
significant  numerical  constants  are  those  that  occur  in  the  heuristic 
expressions  for  the  quanities  a,  3,  and  y. 


TABLE  2.6-1 

Convergence  of  Macroscale  Distribution  in  a  Two-Dimensional  Channel 
with  Parabolic  (Laminar)  Velocity  Profile* 


y/b 

X  /b 
n 

n  ■  0 

n  ■  1 

n  •  2 

n  *  3 

n  -  20 

0.0 

CD 

0.775 

0.695 

0.674 

0.664 

0.2 

DO 

0.775 

0.706 

0.693 

0.689 

0.4 

CD 

0.775 

0.731 

0.727 

0.727 

0.6 

CD 

0.775 

0.758 

0.758 

0.758 

0.8 

CD 

0.775 

0.702 

0.781 

0.781 

1.0 

CD 

0.775 

0.802 

0.800 

0.800 

*  2b  - 

channel  width 

y  ■  distance 

from  centerline. 
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2.7  REYNOLDS  STRESSES  AND  GENERALIZED  EDDY  VISCOSITY  CONCEPT 


In  dealing  with  an  incompressible  fluid,  it  is  customary  to  divide 
the  equations  of  motion  through  by  the  fluid  density.  As  a  result,  the 
viscosity,  pressure  and  stress  terms  all  appear  in  these  equations  divided 
by  the  density.  It  is  convenient  to  avoid  writing  out  the  density  explicitly 
each  time.  Therefore  we  define  the  quotient  of  viscosity  divided  by 
density  as  the  kinematic  viscosity,  and  assign  it  the  symbol  v.  Similarly, 
the  quotient  of  pressure  divided  by  density  becomes  kinematic  pressure, 
symbol  p.  Stress  divided  by  density  becomes  kinematic  stress,  symbol  t^. 

Note  that  kinematic  pressure  and  kinematic  stress  have  dimensions  of 
velocity  squared. 

The  kinematic  Reynolds  stress  tensor  may  now  be  written 

Tij  "  Tij  “  ■5TCIj  (2.7-1) 

The  contraction  obtained  by  equating  indices  and  summing  over  the 
repeated  index  defines  a  scalar  invarient,  namely, 

Tij  "  T11  +  t22  +  t33  -  -<uj  +  **!  +  -  -q2  .  -2E  (-2) 

Now  the  Reynolds  stress  tensor  may  be  rewritten  as  the  sum  of  a 
purely  dllatational  tensor  plus  a  purely  distortional  tensor,  that  is, 

'u  * '  +  Tij  ('3) 

Conversely,  the  distortional  part  is  defined  by 

Tij  '  Tij  +  3"2{ij  ('4) 

Of  course,  the  individual  components  of  depend  on  the  orientation 
of  the  coordinate  exes.  However,  the  double  summation 


1T'  T«  -  t'2 
2  ij  ij 
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is  Invar lent  with  respect  to  any  change  in  the  orientation  of  the  reference 
axes  at  any  point,  and  is  therefore  a  true  scalar.  We  may  think  of  t' 
as  representing  a  kind  of  generalized  scalar  amplitude  of  the  tensor 
Of  course,  the  optional  factor  1/2  was  insetted  into  Eq.  (-5)  to 
maintain  consistency  with  the  corresponding  formula  for  the  strain  rate 
lnvarlent  0.  The  physical  justification  is  as  follows.  Consider  a  simple 
state  of  pure  shear  stress  such  that  “  T22  “  T23  "  T3i  "  °*  Let  ***• 
simple  shear  be  the  only  non-zero  stress.  In  this  case  Eq. 

(-5)  gives 

T'-Ti2  (-e; 

Hence  the  lnvarlent  t*  may  be  identified  with  the  magnitude  of  the  simple 
shear  stress  itself.  Eq.  (-5)  represents  an  extension  of  this  concept  to 
the  general  state  of  stress. 

A  somewhat  similar  analysis  may  be  made  of  the  rate  of  strain  tensor 
of  the  mean  flow  field.  Recall  the  definition 


r 


(-7) 


For  an  incompressible  fluid,  the  associated  dllatlonal  tensor  vanishes, 
that  is, 


,3V  3V  9C, 

rit-ru  +  r22  +  r33-2b^+-s;+ig 

Consequently  the  tensor  is  itself  wholly  dlstortlonal.  An  amplitude 
scalar  (2  may  be  defined  through  the  invariant  relation 


(-8) 


1 

2 


r  r 
U  ij 


(-9) 
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Next  we  wish  to  formulate  an  appropriate  way  to  express  the  general 
relation  that  may  exist  between  the  distortior.al  stress  tensor  and  the 
dlstortional  strain  rate  tensor  .  In  this  connection,  it  is  Instructive 
to  consider  as  a  preliminary  model,  the  relation  which  is  known  to  exist 
between  the  ensemble  averaged  dlstortional  viscous  stress  (t^j)v  and  the 
dlstortional  strain  rates  of  the  mean  flow.  This  is  given  by 

(TI?v  '  vrii  <-10> 

where  v  is  the  ordinary  molecular  kinematic  viscosity. 

This  model  suggests  that  we  attempt  to  express  the  dlstortional 


Reynolds  stresses  in  the  form 

t'  -  el\  . 
ij  13 


(-11) 


where  v  is  now  replaced  by  e,  the  so-called  eddy  kinematic  viscosity. 
However,  this  form  of  equation  implies  that 


.1  Ti  _i  _i  T«  Ti 

T11  a  t22  _  t33  _  t23  _  T31  _  T12 

ril  r22  r33  r23  r31  F12 


(-12) 


Unfortunately,  there  is  no  a-prlorl  reason  to  believe  that  the  tensors 
rjj  and  will  necessarily  satisfy  the  five  similarity  constraints 
expressed  by  Eqs.  (-12).  Ve  therefore  require  a  more  general  relationship, 
one  which  remains  applicable  even  if  and  are  both  specified 
arbitratlly. 

One  method  which  is  sometimes  suggested  for  dealing  with  this 
requirement  is  to  postulate  a  relation  of  the  form 


T  *  ■  r  r 

ij  eijkl*  kl 


(-13) 


Now  the  eddy  viscosity  becomes  a  fourth  order  tensor.  Since  and 
each  have  six  distlnce  components,  the  tensor  has  thirty~slx 
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components  1  By  employing  tensor  arguments,  the  number  of  independent 
eddy  viscosity  components  can  actually  be  reduced  to  ten.  Furthermore,  if 
the  tensors  and  ore  arbitrarily  specified,  Eqs.  (-13)  can  in  fact: 
he  reduced  to  an  equivalent  set  of  five  Independent  equations  in  ten 
unknowns,  however,  these  equations  do  not  constitute  a  determinate  set. 
Clearly,  this  type  of  formulation  is  both  inadequate  and  overly  complex 
as  veil . 


A  somewhat  simpler  tensor  relation  has  been  suggested  by  Dirt, 
Ref.  (4).  In  our  notation  it  becomes 


2l€ikrkj  +  CjkrkJ 


(-14) 


While  this  is  an  Improvement  over  Eq.  (-13),  it  is  not  free  of  difficulties. 
If  is  taken  as  an  unsymmetrlcal  tensor,  Fq.  (-14)  falls  to  produce  a 
determinate  solution  for  when  and  are  specified.  On  the  other 
hand,  if  is  taker,  as  symmetrical,  the  relation  (-14)  loses  generality; 
there  are  then  certain  theoretical  combinations  of  and  which  are 
not  reconcilable  with  (-14)  no  matter  how  the  are  chosen. 

The  foregoing  difficulties  can  be  overcome  by  adopting  the  following 
relation  between  and  1*^,  namely, 

Tij  "  ctrij  +  nfij1  . 


It  will  be  shown  that  the  quantities  c  and  f^  suffice  to  define  the 
relation  between  distortlonal  stress  and  strain  rate  in  a  fully  general 
way.  If  and  are  specified  arbitrarily,  a  corresponding  unique 
solution  can  be  found  for  c  and  f^.  Conversely,  If  c,  f^,  and  are 
specified,  then  t'^  Is  uniquely  determined. 

Solving  Eq.  (-15)  for  f ^  gives 


(-17) 


m  ID*)  wgjewyBnfMrmi1  HWJ  *  » 


It  is  immediately  apparent  from  this  that 


and 


fiJ  "  V 


ii 


fll  +  f 22  +  f33  “  0 


(-n) 


However,  as  it  stands,  Eq.  (-16)  does  not  uniquely  fix  f^  because 
it  contains  the  as  yet  undetermined  parameter  c.  It  will  now  be  shown 
that  for  any  arbitrarily  prescribed  and  r^,  a  corresponding  "preferred" 
value  of  c  can  be  uniquely  defined.  For  this  purpose  let 


t"  ■  efif  ■  t1  —  cT 
ij  ij  ij  ij 

Here  represents  the  discrepancy  between  and  eT^.  We  wish  to 
choose  c  so  as  to  minimize  this  discrepancy  in  some  appropriate  sense. 
Hence  we  form  the  lnvarlent  of  the  discrepancy  namely, 


t"  t"  -  2t 

i^  IJ 


•.2 


lTij  "  erijltTij  ‘  erij3 


-  -  2eFJ  .tl .  +  efcI\T 


IJ  U 


ij  U 


IJ  IJ 


or 


t"2  ■  x'2  -  cT  t'  +  e2n2 
ij  ij 

2 

Now  to  minimize  the  lnvarlent  t"  of  the  discrepancy,  we  set 

,ii2  _ 


3t' 


3e 


-  o  -  r14/ij  +  2en 


whereupon 


(-19) 


(-20) 


(-21) 


r  t  '  r  t  '  r  t 

.  ,  _*  1.1J  .  ij.  IJ 

2  r  r  r  r 

2IT  ij  ij  ij  ij 


(-22) 


This  simple  result  is  important.  It  is  the  basic  definition  of 
the  scalar  eddy  viscosity  c  for  the  case  of  general  stress  and  strain  rate. 
Now  multiplying  (-19)  be  and  using  (-22)  we  obtain  the  significant 


consequence  that 


This  Is  the  mathematical  expression  of  the  concept  that  the  tensors  T 


ij 

and  fjj  are  1  orthogonal" .  Hence  our  basic  stress/straln  rate  law,  Eq. 

(-15),  expresses  the  Idea  that  the  Reynolds  stress  tensor  xj^  consists  of 
a  tensor  which  Is  simply  proportional  to  plus  a  tensor  efif^ 

which  Is  orthogonal  to 

The  tensor  f ^  has  nine  possible  components.  The  symmetry  constraint 

(-17)  reduces  this  to  six  distinct  components.  The  non-divergence  constraint 

(-18)  reduces  it  further  to  five  independent  components.  Tlnally,  the 

orthogonality  constraint  (-23)  reduces  the  degrees  of  freedom  of  f  ^  to  just 

four.  Hence  e  and  f^  together  possess  five  degrees  of  freedom.  This  is 

precisely  the  necessary  and  sufficient  number  to  express  the  general 

relation  between  dlstortional  stress  and  mean  flow  strain  rate.  If  we  add 
2 

as  a  sixth  degree  of  freedom,  this  then  allows  us  to  define  the  general 
relation  htween  overall  Reynolds  stress  and  mean  flow  strain  rate. 

Since  the  foregoing  equations  are  fure  tensor  relations,  they  remain 
valid  for  all  possible  orientations  of  the  coordinate  axes.  It  is  particul¬ 
arly  convenient  to  choose  the  principal  axes  of  strain  rate  in  order  to 
bring  out  the  essential  concepts  in  the  simplest  possible  form.  We  use 
asterisks  to  designate  these  axes.  Also,  there  is  no  loss  in  generality 
in  ordering  the  subscripts  so  that  i  *  r^.  Of  course  * 

r*2  "  0.  Moreover,  it  happens  to  be  useful  to  introduce  the  following 
auxiliary  terminology.  Let 

r*i  ■  r\ 0  "‘“(I’1  +  e* 

r*2  -  fi  8ln(0  +  e*>  (-24) 

r53  ■  A 0  ¥ + H 
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The  merit  of  Eqs.  (-24)  lies  in  the  fact  that  they  express  the  three 
quantities  T *^ ,  in  terms  of  just  tvo  parameters,  namely,  ft  and  0* 

Furthermore,  as  the  reader  may  easily  verify,  eqs.  (-24)  satisfy  identically 
the  following  important  constraints, 


r*  +  r*  +  r*  ■  o 
11  22  *33 


2  2  2 
r*r  +  r*r  +  r** 

'll  22  1  33 


-  ft 


The  parameter  0*  is  then  fixed  from  the  relation 


(-25) 


j  3  r* 
v  1  22 


(r!l  "  *$3* 


■  tan  0* 


(-26) 


It  is  eqislly  useful  to  adopt  a  similar  scheme  for  the  f .  Thus 
f •j=a|  g*  sin|~  +  <|i* 


whereupon 


f 22  "77  8*  8ln^°  +  **) 
f33  m  A  R*  8ln(~  T~  +  **] 


f*  +  f*  +  f*  .  o 
11  r22  *33 


(-27) 


(-28) 


2*>? 

f*"  +  f*-  + 

11  *22  *33 


m  p* 


and 


H  lh 

(fii  -  {y 


■  tan  ip* 


(-29) 
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Now  the  orthogonality  constraint,  er.  (-23),  reduces  in  this  case  to 


the  form 


r*  f*  +  r*  f*  +  r*  f*  m  n 
1 llrll  1 22r22  1 33r33 


(-30) 


Upon  substituting  Eos.  (-24)  and  (-27)  into  (-30)  and  reducing,  we 
obtain  the  beautifully  simple  result 


¥  -  6*  ±  |  (-31) 

This  too  expresses  a  kind  of  orthogonality  between  the  and  the  f^. 

With  ip*  now  known,  the  components  f ,  f^»  *§3  are  all  fixed  as 
soon  as  g*  is  specified.  The  off-diagonal  terms  f^,  f^,  and  f*^  may  be 
specified  individually.  The  final  stress/straln  rate  equations  for  these 
particular  axes  now  taken  the  form 
“7  12 

Tfi  ■  '“I  ■  '  3  11  +  tlrh  +  “'h1 
t22  -  -»2  -  -  f  I2  ♦  «irj2  ♦  t-f|2) 


r#r 

•33 


■  -  i  "2  +  ‘r53  +  "fei 


(-32) 


t23  *  ~U2U3  "  0  +  eI°  + 


*  ■  -u^u^  ■  0  +  c[0  +  fif*,,) 


T31 


31 


t12  ‘  '“lu2  ‘  0  +  E<°  +  0f5l] 

Mote  that  the  principal  axes  of  stress  and  strain  will  coincide  if, 
and  only  if,  f*^  -  fjj^  -  fjj  "  0. 

The  further  reduction  of  these  eauations  for  the  important  special 
case  of  a  plane  mean  flow  is  interesting.  Now  the  axes  x^  and  x^  are  in  the 
plane  of  the  flow,  and  X2  is  normal  to  this  plane.  Hence  we  have 
0*  •  0 


ril  -  HI 


(-33) 


82 


Consequently  from  Sq.  (-31)  we  find  that 


.  j.  .11 

V*  *  ±  2 


(-34) 


There  la  some  choice  concerning  the  manner  of  handling  algebraic 
signs.  We  choose  to  drop  the  negative  sign  in  (-34)  and  instead,  to  allow 
g*  In  the  following  expressions  to  tahe  on  positive  or  negative  values  as 
required.  Then 


ftl  “  -  7?  6* 


fh-o 


f*  m  + 

a  n  » 


22  " 


f^  •  arbitrary 


(-35) 


f*3  -  g< 


f*  .  o 

l2 


Of  course,  f*3  ■  f*2  -  0  by  reason  of  symmetry,  hence  p*  and  f^ 
remain  as  the  two  arbitrary  degrees  of  freedom  which  define  the  coefficients 


The  final  atreas/strain  rate  equations  become 


Ttfc  ■ 
T11 


-  3  q2  +  -  -5  R*1 


r  ft  m 
T22 


-  \  q2  +  ec[0  +  -|  g*J 


-  I  o2  +  ef![-l  -  g*J 


(-36) 


t31  "  ”U3U1  “  U  +  efif31 


r#f  m  rtfc  m  0 

t23  t31 


Hence  the  overall  strrss/strain  rate  relation  is  now  fully  specified 
1  2 

by  the  four  quantities  q  ,  e,  g*,  and  f*^. 
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In  Che  case  of  axi-symmetric  mean  flows,  the  results  are  similar  to 
those  for  plane  flow,  except  that  G*  Is  not  In  genera]  zero;  therefore  ip* 


must  be  found  from  Eq.  (-31)  and  the  f*j  from  Eq.  (-27).  However,  we  still 


have  •  f*^  «  0  for  this  case. 


1  2 


The  basrc  strcss/straln  parameters  are  now  -j  q  ,  e,  g*,  f^,  and  \p*. 
Note  from  Eqs.  (-20  and  (-31)  that  ip*1  Is  a  known  function  of  the  strain 
tensor,  but  that  the  other  four  parameters  may  be  arbitrary. 

For  purposes  of  approximation  in  connection  with  our  heuristic 


turbulence  model,  we  have  decided  to  make  the  simplifying  assumption  f 


ij 


0. 


Actually  this  Is  not  as  drastic  as  it  seems.  In  the  first  place,  as  we 
have  seen,  f^  embodies  at  most  only  four  degrees  of  freedom,  not  six  or 
nine.  Secondly,  we  are  mainly  concerned  with  plane  or  axisymmetrlc  mean 
flows,  for  this  f^,  has  only  two  degrees  of  freedom,  namely,  g*  and  f^. 

Thus  this  assumption  amounts  to  neglecting  but  two  degress  of  freedom. 

Secondly,  an  analysis  of  thu  momentum  transfer  Involved  in  generating 
the  Reynolds  stresses  suggests  that  the  term  accounts  for  most  of  the 
effect,  and  that  the  deviation  efifjj  should  Indeed  he  small.  While  the 
experimental  data  which  Is  available  Is  not  decisive  on  this  point,  it  does 
tend  to  support  the  above  hypothetic.  It  does  seem  that  this  simplification 
may  become  somewhat  Inadequate  In  the  Immediate  vicinity  of  a  fixed  wall, 
but  it  appears  to  be  reasonably  satisfactory  elsewhere  in  the  flow. 

Under  these  assumptions,  the  principal  axes  of  stress  and  strain 
coincide,  and  the  equations  for  these  particular  axes  reduce  to  the  simple 
form 


1  _2 
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T11  "  “  ^  q~  +  efi 
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The  foregoing  stress/strain  rate  law,  despite  its  apparent  simplicity, 
can  be  useful  for  calculating  the  Reynolds  stresses  which  correspond  to 
any  specified  mean  flow  strain  rates  only  if  the  value  of  e  can  somehow  he 
estimated  with  acceptable  accuracy.  Whether  this  can  be  done  is  an 
empirical  question,  capable  of  being  answered  by  suitable  experimental 
information.  The  experimental  data  so  far  available  is  not  yet  adequate 
to  answer  this  need  fully.  Further  experiment  and  analysis  is  needed. 
Tentatively,  however,  we  are  justified  in  postulating  that  c  is  some  function 
of  the  characteristic  local  turbulence  parameters  q,  X,  and  t  as  prevlsouly 
defined,  perhaps  also  of  molecular  kinematic  viscosity  v,  distance  y  from 
the  nearest  wall,  and  possibly  of  some  as  yet  unknown  additional  parameters 
as  well.  This  can  be  written  symbolically  in  the  form 

e  -  e(q,  X,  t,  v,  . . )  (-38) 

Actually,  we  would  be  on  slightly  more  fundamental  ground  to  choose 
X*  and  t*  in  place  of  X  and  t.  Unfortunately,  X*  and  i*  are  not  usually 
known,  whereas  X  and  t  have  been  related  in  a  specified  manner  to  the  mean 
flow.  Consequently,  for  pragmatic  reasons  we  replace  X*  and  t*  by  X  and  t 
in  the  above  expression. 

The  foregoing  empirical  law  which  governs  the  eddy  viscosity  e  can  be 
greatly  simplified  by  dimeaslonal  analysis.  We  recall  that  in  mechanics 
generally,  the  number  of  fundamental  dimesnions  is  three.  These  are  often 
taken  as  mass,  length  and  time.  However,  in  the  case  of  incompressible 
flow,  dlvldlon  of  the  equations  of  motion  by  density  has  the  effect  of 
eliminating  mass  from  the  problem.  Our  final  equations  are  therefore  all 
klnematlcal  in  character,  with  length  and  time  as  the  only  fundamental 
dimensions.  It  follows  therefore  that  any  two  of  the  above  parameters  may 
be  chosen  as  the  necessary  and  sufficient  basis  for  the  non-dimensional izing 
process. 
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From  a  purely  logical  standpoint,  it  would  seem  natural  and  convenient 
to  use  the  basic  local  scales  of  length  and  time,  A  and  t,  for  this  purpose. 
However,  there  are  grounds  for  believing  that  q  and  A  play  roles  which  are 
physically  more  important  than  that  played  by  t.  This  seems  to  be  especially 
true  for  mean  flows  which  are  steady,  or  which  vary  only  slowly.  Conse¬ 
quently,  we  choose  q  and  A  and  our  non-dimenslonalizlng  parameters. 

As  a  result  of  the  foregoing  considerations,  the  empirical  eddy 
viscosity  law  may  now  be  rewritten  in  the  dimensionless  form. 


fiJ  .  ,  .  .ffll  fli  X  I 

IqAj  -  °[a  *  v  *  A  *  *  *  * J 


(“39) 


Thus,  our  problem  has  been  simplified  to  that  of  finding  the  dimen¬ 
sionless  coefficient  a  which  presumably  varies  as  a  function  of  the  dimen¬ 


sionless  parameters 


fail 


and 


^1  and  [? 


,  and  possibly  also  of  one  or  two 


other  dimensionless  variables  as  well. 

In  practice,  the  usefulness  of  the  above  formulation  hinges  on  whether 
the  dimensionless  coefficient  a  does  in  fact  behave  in  a  reasonably  simple 
and  stable  manner.  Our  present  indications  are  that  this  is  indeed  so. 

Owing  to  limitations,  of  time  and  data,  we  have  not  fully  established  the 
actual  dependence,  if  any,  of  a  on  the  quantities  |^j  and  .  We  have 
four,  heuristlcally,  however,  that  satisfactory  overall  results  can 
apparently  be  found  by  assuming  that  a  depends  only  on  the  wall  distance 
parameter  (y/A) .  The  following  expression  seems  to  give  reasonably  good 
curve  fits,  namely, 

*12, 


-  0.065  £ 


1  +  e 


I- 1 
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(“40) 


Of  course,  this  expression  is  tentative.  More  data  is  needed  on  this 
point  in  any  ca6e. 
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Naturally,  this  particular  expression  for  a  Is  applicable  only  so 
long  as  we  specifically  retain  our  present  definition  of  A.  Any  modifica¬ 
tion  of  the  definition  of  A  clearly  calls  for  a  corresponding  adjustment 
in  a.  Note  however  that  the  dimensionless  quantity 


should  remain  substantially  invarlent  with  respect  to  such  details  * 

However,  even  if  a  be  accurately  known,  and  X  as  well,  the  actual 

eddy  viscosity  e  cannot  be  evaluated  until  q  also  is  established.  Bearing 

1  2 

in  mind  that  the  turbulent  energy  is  E  =.  a  ,  it  is  apparent  that  the 
determination  of  q  must  involve  an  energy  analysis  of  the  turbulence. 
Physically,  this  is  fundamental.  It  ties  together  the  turbulence,  the 
eddy  viscosity  and  the  mean  flow  into  a  consistent  overall  picture.  It 
also  provides  the  additional  equations  which  are  necessary  to  form  a 
determinate  set. 

2.8  HEURISTIC  ENERGY  ANALYSIS  OF  THE  TURBULENT  FIELD 


Starting  with  the  Navier  Stokes  equations  of  motion,  it  1b  possible 
to  develop  a  corresponding  exact  energy  equation,  as  shown  in  detail  in 
Appendix  A.  The  following  result  is  obtained,  namely, 


2  ±j  i j 


II  III 
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V  (2.8-1) 
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This  equation  governs  the  evolution  over  time  and  distribution  over 
space  of  the  ensemble  mean  turbulent  kinetic  energy,  that  is,  of 

E  '  1  1  -  2  °juj  '  2  <“l  +  u2  +  “3>  <-» 

The  five  terms  on  the  right  side  of  Eq.  (-1)  have  the  following 
physical  significance: 

Term  I  -  The  rate  of  energy  input  into  turbulence  as  work  done  by  the  mean 
flow  against  the  resistance  of  the  Reynolds  stresses. 

Term  II -The  rate  of  irreversible  energy  dissipation  into  heat.  Note  that 
the  quantity  Is  positive  definite.  This  quantity  accounts 

for  the  local  entropy  production. 

Term  III-The  r s\.e  of  convective  transport  of  turbulent  energy  F  by  the  mean 
flow  U . . 

Term  IV  -The  rate  of  net  transport  by  turbulent  diffusion  of  the  turbulent 
pressure  plus  kinetic  energy  fluctuations. 

Term  V  -  The  rate  of  net  transport  of  energy  by  molecular  diffusion.  Note 
that  the  Reynolds  stresses,  as  well  as  the  energy,  appear  in  this 
term. 

While  Eq.  (-1)  is  theoretically  exact,  it  contains  several  terms 
which  cannot  he  evaluated  exactly  without  detailed  knowledge  of  the 
turbulent  fluctuations.  Such  detailed  information  is  not  practically 
attainable.  It  is  therefore  necessary  to  approximate  these  terms  by  suit¬ 
able  semi-empirical  expressions.  These  expressions  will  necessarily  contain 
empirical  coefficients  whose  values  can  in  principle  be  chosen  so  as  to 
optimize  agreement  with  available  experimental  data. 

We  have  shown  in  the  preceding  section  that  the  Reynolds  stresses 
can  always  be  expressed  with  complete  exactness  and  generality  in  the  form 
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(-3) 
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where 


e  • 


Vaj 
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and 


(-4) 


(-5) 


Consequently,  term  I,  the  work  Inrut  expression  in  the  energy  equation 
reduces  exactly  to  the  result 


2Tijrij  S  cft 


(-0 


however,  as  we  have  seen  earlier,  c  must  here  be  approximated  by  an 
expression  of  the  form 

e  ■  aXq  (-7) 

where  a  is  presumably  a  well  behaved  empirical  function,  the  expression 
for  which  has  been  given  earlier. 

Next  consider  the  dissipation  of  energy  into  heat,  as  given  by  term  II 
of  the  basic  energy  equation.  Dimensional  considerations  suggest  that 
this  term  can  be  expressed  to  good  advantage  in  the  form 


2VYijYij 


(-P) 


This  expression  amounts  to  a  definition  of  the  so-called  dissipation  length 
Ap  which  is  a  length  characteristic  of  the  microscale  of  turbulence.  This 
form  of  expression  would  be  very  convenient  if  the  dissipation  length  A^ 
happened  to  be  related  to  the  macroscopic  scale  A  in  some  fairly  simple  and 
invar lent  way.  Unfortunately,  this  does  not  seem  to  be  the  case. 
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However,  experimental  studies  such  of  those  by  kaufer,  Pef.  (6), 
clearly  suggest  that  at  sufficiently  high  Reynolds  numbers,  the  heat 
dissipation  effects  tend  to  become  independent  of  Reynolds  number!  This 
Implies  that,  contrary  to  Eq.  (-8),  the  optimal  form  of  expression  for  Eu 
is  one  which  does  not  contain  viscosity  v  explicitly.  This  will  be  the 
case  in  fact  if  we  asstime,  for  example,  that 


where  61  is  now  hopefully  some  well  behaved  empirical  function.  For 
combining  these  last  two  results  gives 


h 


(-10) 


which  is  not  explicitly  dependent  on  v.  This  form  is  similar  to  that 
suggested  by  Prandtl  in  Fef.  (2). 

The  fact  that  the  viscous  dissipation  occurs  in  a  form  which  does  not 
explicitly  contain  the  viscosity  seems  rather  paradoxical  at  first  sight, 
and  calls  for  some  comment.  It  can  be  given  a  reasonable  interpretation  by 
considering  turbulence  energy  processes  from  the  spectral  point  of  view. 

It  is  well  known  that  the  work  input,  term  I,  occurs  mainly  at  the  long 
wave  length  end  of  the  spectrum,  and  that  the  dissipation  to  heat,  term  II, 
occurs  mainly  at  the  short  wave  length  end.  At  high  Reynolds  numbers, 
there  is  a  broad  range  of  wavelengths  in  between,  the  so-called  inertial 
range,  for  which  both  work  input  and  heat  dissipation  are  negligible.  Each 
wave  length  within  this  inertial  range  simply  receives  energy  from  the 
longer  wave  length  components  and  transmits  this  energy  to  the  shorter 
wave  length  components.  It  would  appear  that  the  rate  of  transmission  of 
energy  thru  the  spectrum  is  largely  controlled  by  this  process  in  the 
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Inertial  range,  rather  than  by  the  viscous  dissipation  process  Itself. 
Apparently  the  viscous  dissipation  rate  easily  adjusts  Itself  as  nay  be 
required  to  dissipate  all  of  the  energy  coning  through  the  inertial  range. 
Since  the  viscous  process  Is  not  rate  controlling,  the  effect  of  viscosity 
tends  to  drop  out  of  the  experlnental  picture. 

These  experimental  facts  and  theoretical  considerations  suggest 
that  the  heat  dissipation  should  be  analyzed  by  the  same  general  method 
previously  employed  In  connection  with  the  eddy  viscosity.  For  example, 
we  might  assume  tentatively  that 


Ey  -  Ey(q,  X,  t,  v,  y,....)  (-11) 

Then  non-dlmcnslonallzlng  on  the  basis  of  q  and  X  as  common  parameters,  we 
obtain 


(-12) 


which  is  seen  to  be  equivalent  to  Eq.  (-10) 

This  result  suggests  that  £'  is  a  well  behaved  universal  function 
which  can  be  established  empirically.  However  sound  In  principle,  this 
idea  is  rather  difficult  to  apply  at  the  present  time  owing  to  the  paucity 
of  suitable  experimental  data.  It  does  seem  plausible  that  for  mean  flows 
which  are  not  too  strongly  unsteady,  the  effect  on  B'  of  the  parameter 

might  prove  to  be  negligible.  On  the  other  hand,  efforts  to  establish 
6*  as  a  function  of  alone  have  not  as  yet  proved  conclusive  one  way 

or  the  other.  This  Is  due  In  part  to  limitations  in  the  data  and  In  part 
to  limitations  In  time. 

Instead  of  following  the  above  line  of  Investigation  to  the  end,  we 
discovered  the  following  more  promising  modification.  Numerical  experimenta¬ 
tion  suggested  that  It  might  provide  a  more  satisfactory  alternative. 
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Recall  that  the  quantities  A  and  t  In  our  model  are  themselves  computed 
from  the  Integrals  I  and  J.  Hence  It  amounts  merely  to  a  rearrangement  of 
form  to  rewrite  F.q.  (-11)  as  follows: 

"  Hr ( q »  I,  J*  v,  . . )  (—13) 

Nov,  non-dimensional Izlng  on  the  basis  of  q  and  J  as  common  parameters,  we 
obtain 


«,7j)1/3 


(-14) 


Actually,  the  significance  of  the  dimensionless  group  on  the  left  was 
ertabllshed  by  certain  somewhat  speculative  physical  arguments  which  will 
not  be  repeated  here.  (See  ref.  11)  It  was  found,  however,  on  the  basis  of 
numerical  Investigation,  that  this  particular  dissipation  parameter  0  seems 
to  remain  more  nearly  constant,  and  to  agree  better  with  physical  data 
than  docs  the  parameter  S'  resulting  from  the  earlier  analysis.  Hence 
this  Is  the  form  that  was  retained  for  our  present  heuristic  turbulence 
model. 


Here  0  Is  a  slowly  varying  dimensionless  function.  The  dependence 

1  -^1/3 


of  0  on  the  quantities 


(qJ) 


2/3  and  v 


J 


was  not  actually  Investigated, 


due  to  time  limitations.  However,  It  was  found  heurlstically  that  0  seems 

fil 


to  depend  only  on  the  relative  distance  ^-jto  the  nearest  wall,  If  any, 
and  that  far  from  any  walls  0  becomes  essentially  constant. 

The  followli  g  empirical  equation  seems  to  produce  reasonably  good 
curve  f Its ,  namely. 


e 
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(-15) 
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Some  data  which  tends  to  confirm  this  formulation  is  given  in 
Appendix  C.  In  any  case,  It  Is  clear  that  the  energy  dissipation  to  heat 
needs  further  investigation. 

Next  consider  term  IV  of  the  basic  energy  equation.  It  la  useful 
for  purposes  of  discussion 


P  + 


V 


M  2 


(-16) 


The  second  of  these  terms  clearly  represents  the  net  transport  of  the 

(UJUll 

fluctuating  turbulent  kinetic  energy  by  the  action  of  the  turbulent 

fVll 

velocity  fluctuations  u^.  Clearly  the  turbulent  kinetic  energy  I-  ■"-*1  is 
both  a  scalar  and  an  extensive  property.  It  may  be  regarded  as  something 
which  is  physically  transported  along  with  the  element  of  fluid  mass.  In 
this  respect  it  is  quite  analogous  to  other  extensive  properties,  such  as 
salinity,  for  example.  It  is  customary  and  appropriate  to  represent  such 
net  turbulent  transport  by  an  appropriate  empirical  diffusion  coefficient. 
Hence  we  write 


(-17) 


On  the  other  hand  the  pressure  velocity  correlation  term  u^p  Is  in 
a  somewhat  different  category  owing  to  the  fact  that  pressure  p  is  an 
intensive  property.  It  is  not  an  extensive  property  which  can  be  regarded 
as  being  transported  along  with  the  fluid  mass.  We  note  however  that  p  is 
a  dependent  variable.  This  means  that  whenever  the  spatial  distribution 
of  the  velocity  fluctuations  u^  is  specified,  the  corresponding  spatial 
distribution  of  p  is  likewise  fixed.  Hence,  looking  at  the  matter  statisti¬ 
cally,  we  can  say  that  p  is  correlated  in  some  fashion  with  u^.  In  fact, 
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taking  Into  account  dimensional  considerations,  we  can  improve  this  state* 

fVil 

ment  and  say  that  p  is  statistically  correlated  with  - 1  .  For  example, 
in  the  limiting  case  of  an  inviscid  fluid,  Bernoulli’s  equation  informs 
us  that  regions  of  higher  than  average  kinetic  energy  will  tend  to  be 
regions  of  lower  than  average  pressure,  and  vice  versa.  Consequently  we 
are  justified  in  writing  the  pressure  term  in  the  analogous  form 


Now  the  overall  effect  becomes 

ul(,,+-2-1)  ■■<C"  +  e',(%}  <-W> 

normally  c'  is  a  positive  quantity,  however,  because  of  the  generally 
negatfye  correlation  befveen  pressure  and  kinetic  energy,  we  expect  e"  to 
be  usually  negative  and  smaller  than  c'.  Bence  (e"  +  e')  should  be 
positive  but  smaller  than  c'  alone.  Moreover,  e"  and  r.'  are  both  of  the 
same  dimension  as  the  kinematic  eddy  viscosity  c.  Also  e'  can  be  expected 
to  be  roughly  of  the  same  order  of  magnitude  as  c.  For  these  reasons,  it 
is  advantageous  to  write. 


(c'  +  c')  -  ye  (-20) 

where  y  is  not  a  dimensionless  coefficient  roughly  of  order  of  magnitude 
uni ty . 


Hence  finally  we  write 


r 

u.  p  + 
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-ye 


3E 


3x , 


(-21) 


Various  numerical  experiments  suggest  that  the  dimensionless  empirical 
coefficient  y,  like  the  coefficients  a  and  f,  depends  only  on  the  distance 
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y  to  the  nearest  fixed  boundary,  If  any.  The  followinp.  expression  has  been 

found  to  give  reasonably  good  results. 

,2 


Y  ■  1.4  -  0.4e 


-fr  -  >1 


(-22) 


Turning  next  to  tern  V  of  tht  basic  energy  equation,  the  molecular 
diffusion  effect,  It  nay  be  said  that  this  will  normally  be  very  much 
smaller  than  the  other  terms  In  the  equation.  It  may  safely  be  neglected 
everywhere  except  In  the  very  thin  laminar  sublayer  region  along  any  solid 
boundary.  We  shall  Ignore  It  In  our  mode] . 

Finally  we  combine  the  foregoing  steps  to  obtain  an  overall  heuristic 
energy  equation  as  follows.  Let 


2 

q*  -  2E 


Then  the  energy  equation  becomes 
(ft)  “  cB2-e(2E)7/6Jl 


(-23) 


(-24) 


where  c,  6,  end  a  are  fixed  by  the  following  empirical  expressions,  namely, 

l  2 


0.065  jl- 
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12.  Batchelor,  G.  K.,  "Theory  of  Homogeneous  Turbulence",  Cambridge 
University  Press,  1953. 
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13.  Pal,  S.  1.,  "Viscous  Flow  Theory:,  Van  Nostrand,  Frlnceton,  1956. 

14.  Hlnxe,  J.  0.,  "Turbulence  -  An  Introduction  to  Its  Mechanism  and 
Theory",  McGrav  Hill,  1959. 

For  specific  accounts  of  more  recent  research,  perhaps  the  two  most 
useful  journals  are  the  "Journal  of  Fluid  Mechanics"  and  the  "Physics  of 
Fluids." 
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APPENDICES 


APPENDIX  A 


Derivation  of  Exact  Energy  Equation  for  Turbulence 


Equation  of  Notion 

Let  prime  marks  on  a  symbol  be  used  to  denote  the  sum  of  the  mean 
value  of  that  quantity  plus  the  turbulent  perturbation.  That  is,  let 


j 


ui  + 


P'  -  P  +  p 

Also  let 


(A-l) 


(-2) 


denote  viscous  shear  stresses. 

The  equations  of  continuity  and  motion  now  become 


(-3) 


(-4) 


The  corresponding  energy  eouation  is  found  by  multiplying  the  equation 
of  motion  by  Uj.  Thus 


Now  we  substitute  Eqs.  (-1)  and  (-2)  into  (-3),  (-4)  and  (-5),  then 
expand,  average  and  simplify.  All  terms  linear  in  perturbation  quantities 
vanish  in  the  averaging.  Therefore,  the  averaped  equations  of  continuity, 
motion  and  energy  become,  respectively, 
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0 


(-6) 


£ i 

U*i 


Also,  multiplying  the  averaged  equation  of  motion 
plves  the  energy  equation  for  the  mean  flow,  namely. 


(-7)  through  by  l'j 


-  l\ 


3 

3*. 


<LiU3  +  uiu3> 


Now  subtracting  (-10)  from  (-9)  gives  the  energy  equation  for  the 
turbulence  Itself,  that  Is 


. .  eld .  »_r. 

i]  .  3  eaa 

3tl  2  ) 

j13«1  J  3*,L 

1  2  J  UiUj  [3* 

4J  a*il  2  J 

It  Is  convenient  to  symbolize  the  mean  turbulent  kinetic  energy,  the 
Instantaneous  kinetic  energy,  and  the  Reynolds  stresses,  respectively,  by 

( 
( 

( 

Therefore  the  energy  equation  becomes 


3^(U1E)  +  Tij 


( 
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-7) 

-8) 

-9) 

-10) 

-11) 

-12) 

-13) 

-14) 


The  viscous  tern  can  now  be  developed  further  as  follows.  Consider 


the  identity 


where 


il  3  rvi 

l)  *  rx7<Wi>  -  *u(kJ| 

(-15) 

the  viscous  stress  perturbation  nay  be  written 

*1J  "  ' 

+  VYu 

(-If) 

r3u ,  Su.^ 

YU" 

_ i  + _ i| 

l3*i  3*J 

(-17) 

and 


Therefore 


-  -utp  +  -jjji] 


S3--Sij 


0  +  v  Yd 

Substituting  (-19)  and  (-20)  into  (-15)  gives 


Sf) '  +  v(%  ^fl]  -*»u^ 

Substituting  this  back  into  (-14)  and  rearranging,  we  find 

bd  .  ,  P]  .  vTpi  -  -JL(l.  fj 
UtJ  TiJia*1J  vriJi3*1J  ax1lliEJ 


3 

3x, 


(-18) 


(-19) 


(-20) 


(-21) 


(-22) 


101 


Two  further  simplifications  are  advantageous.  Recall  that 


-  -  CHj  +  ~k\ 

ij  axj) 


(-23) 


Noting  further  that  and  are  both  synnetrlcal  tensors,  ve  write 


or 


.  T  G±  +  -  2t  [!!i| 

ij  "ij^xj  axjj  "ijUxJ 


&1  -  I,  r 

ij^xj  2Tij‘ij 


By  an  exactly  analogous  irgument  ve  find  that 


nsi  -w- 

*  ijYij 


Consequently,  the  turbulence  energy  equation  becomes  finally. 


(ft) "  rijrij  ■  r^ITTj "  ^(uiE) "  ui(p  +  -V)] 


“  fUL 

_  Hull 

.  '3xi 

^  J 

(-24) 


(-25) 


(-26) 


(-27) 


It  is  Important  to  notice  that  this  equation  is  complete  and  exact. 
Each  term  has  a  definite  physical  significance  as  discussed  in  the  main 
body  of  the  report.  Th*  method  of  evaluating  certain  of  the  terms  heurls- 
tlcally  is  also  discussed  in  the  main  report. 
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APPENDIX  B 


Local  Length  Scale  for  Axi-Symmetric  Hear  Flowq 


Consider  the  following  sequence  of  relations  pertaining  to  the 
heuristic  model,  where  the  various  integrals  are  summed  over  the  entire  flow 


field. 


_  AS* Ax 
xi 

W  ■  e  ■  Weighting  function 

H  ■  /Wdv'  ■  Normalizing  integral 
W 

w  ■  —  ■  Normalized  weighting  function 

2  4 

I  **  /wfi  dv'  *  First  characteristic  Integral 
2  2 

J  ■  /w(fiD')  dv'  -  Second  characteristic  integral 
2  1^ 

X  ■  — r-  ■  Local  scale  parameter 


(B-l) 

(-2) 

(-3) 

(-4) 

(-5) 

(-6) 


For  axi-symmetric  flows,  the  volume  Integrals  (-2),  (-4),  and  (-5) 
reduce  to  surface  integrals  over  the  meridional  plane. 

In  this  case  the  volume  element  may  he  written  in  the  form 
dv'  -  r'dtf'dr'dz' 

Also  in  cartesian  coordinates 


Ax*Ax  -  (x1  -  x)2  +  (y'  -  y)2  +  (z'  -  z)2  (-7) 

We  may  express  the  coordinates  of  the  fixed  point  x  and  of  the 

r 

variable  point  X'  in  terms  of  cylindrical  coordinates  as  follows. 
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X  ■ 

0 

X* 

-  r' 

sin 

e’ 

y  ■ 

r 

y’ 

-  r' 

cos 

e’ 

(-8) 

z  ■ 

z 

z' 

-  z’ 

I'pon  substituting  (-8)  into  (-7)  and  rearranging,  we  can  write 

K  -  C2  +  2rr’(l  -  cos  6')  (-9) 


where 

C2  *  <r’  -  r)2  +  (*'  -  z)2  (-10) 

2 

Now  let  us  consider  the  typical  Integral  I  as  defined  by  Eq.  (-4). 

By  making  use  of  the  foregoing  relations,  this  Integral  csn  he  reduced  to 
the  required  form.  In  addition,  it  is  convenient  to  utilize  the  follovrlnp 
auxiliary  variable,  namely, 


2 


(-11) 


Eq.  (-4)  is  now  expressible  in  the  form 


How  since  f2  la  a  function  of  z'  and  r'  only,  but  independent  of  0', 
the  above  integration  with  respect  to  6'  can  be  made  once  and  for  all. 
This  leads  naturally  to  the  definition  of  the  auxiliary  function 


*(0  -  ^  J  e  "2c2(1  "  ct>8  e  )de'.  (-13) 

o 

The  function  ij«  defined  above  muBt  be  computed  numerically.  The 
factor  — *■  was  Inserted  on  the  right  so  as  to  achieve  the  simplification 
that  at  large  values  of  £  the  function  i*  approaches  the  limit  unity,  that 
Is,  that 
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lim  ip  ■  1 


(-14) 


It  then  happens  that  for  small  values  of  t 

£im  iji  ■  2/Tr  c  (-15) 

;-*o 


The  function  i|<  is  shown  in  Fig.  R-] .  Computed  values  are  listed  in 
Table  B-l. 

2 

Once  the  function  is  known,  the  integral  1  may  be  rewritten  as 
follows. 


Similarly 


(-17) 


2 

The  scale  parameter  A  is  now  found  from  (-f ) .  Notice  that  in 
2  1  /tT 

computing  A  from  (-6)  the  factors  ^  J—  in  (-16)  and  (-17)  cancel  from  the 
result. 

Eqs.  (-16)  and  (-17)  illustrate  how  the  characteristic  integrals 
reduce  from  volume  to  surface  integrals  for  the  axi- symmetric  case.  Eqs. 
(-11)  and  (-13)  define  i|>  which  is  a  function  characteristic  of  axi-symmetric 
mean  flows. 
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TABLE  B-l 


COMPUTED  values  of  THE  y  functiok 
(Refer  to  Eqs.  B-ll,  b-13) 


0.0 


0.1 


0.2 

0.3 

0.4 

0.5 

0.6 

0.7 

0.8 

0.9 

1.0 

1.2 

1.4 

1.6 

1.8 

2.0 


2.5 


3.0 


4.0 


5.0 


0.0000 

0.3475 

0.6555 

0.8955 

1.0562 

1.1433 

1.1739 

1.1687 

1.1461 

1.1187 

1.0936 

1.0581 

1.0387 

1.0279 

1.0213 

1.0169 

1.0105 

1.0072 

1.0040 

1.0025 

1.0000 


107 


APPENDIX  C 


Experimental  Data 

Various  useful  reference  Items  of  experimental  data  are  summarised  In 
Tables  C-l  through  C-5  and  In  Figs.  C-l  and  C-2.  These  tabulations  and 
dlagra/ns  are  largely  self-explanatory. 

Table  C-2  is  of  special  Interest  In  that  it  shows  a  comparison 
between  theoretical  values  of  dissipation,  as  computed  from  the  heuristic 
model,  and  the  experimental  measurements  of  Laufer,  Ref.  (10).  The  agree¬ 
ment  Is  satisfactory. 

Laufer  *8  measured  values  are  known  to  be  too  low  and  he  attempted  to 
extlmate  better  values  as  shown. 

Shearing  rates  were  estimated  on  the  basis  of  a  theoretical  mean 
velocity  profile  using  Prandtl's  mixing  length  theory  and  Nlkuradse's 
empirical  curve  for  the  mixing  length  t,  namely, 

l  -  0.1A  -  0.08r2  -  0.06r4  (C-l) 

The  Information  In  Table  C-2  Is  also  displayed  graphically  In  Figs. 

C-l  and  C-2.  There  Is  no  particular  significance  to  the  fact  that  the 
ordinate  in  Fig  C-2  is  ^  rather  than  B  Itself;  this  merely  reflects  the 
circumstance  that  a  slightly  different  notation  was  being  employed  at  the 
time  these  calculations  were  actually  made  than  that  which  is  used  In  this 
report . 

In  view  of  the  uncertainties  In  the  experimental  data  itself,  the 
agreement  between  theory  and  experiment  as  revealed  in  Figs.  C-l  and  C-2 
is  regarded  as  satisfactory. 
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The  remaining  tables  require  no  particular  comment  except,  perhaps, 
to  note  that  the  available  data  Is  meagre  and  not  very  satisfactory.  Note 
that  the  data  given  Is  mostly  one  dimensional,  with  slight  information  In 
the  second  dimension  and  none  at  all  In  the  third.  There  is  clearly  a 
need  for  more  experimental  vcrk. 


TABLE  C-l 


TURBULECT  ENERGY  DATA  FOR  A  PIPE 


Measurements  by  Laufer,  Ref.  (10) 


All  quantities  non-dlmenslonallzed  using  radius  a  ■  1  and  friction 
velocity  v*  ■  1. 


R  - 

e 

500,000 

Re  -  50,000 

r 

2E 

h 

2E 

• 

■to 

Ti/ 

0 

1.65 

2.1 

1.69 

’.15 

— 

0.1 

1.80 

2.2 

1.80 

2.26 

— 

0.2 

2.05 

2.4 

2.06 

2.72 

— 

0.3 

2.45 

3.1 

2.41 

3.28 

1.5 

0.4 

2.90 

3.8 

2.85 

3.96 

2.5 

0.5 

3.50 

4.8 

3.48 

4.76 

3.7 

0.6 

4.20 

6.2 

4.01 

6.23 

5.3 

0.7 

5.00 

8.7 

4.49 

8.10 

7.6 

0.8 

6.00 

13.0 

5.13 

11.3 

11.7 

0.9 

7.40 

22.8 

5.66 

20.9 

22.0 

0.98 

— 

— 

7. 85 (peak)  — 

— 
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TABLE  C-2 


MEASURED  AND  COMPUTED  ENERGY  DISSIPATION 
RATES  IN  A  PIPE  AT  Rft  -  500,000 

Measurements  by  Laufer,  Ref.  (10) 


A1  quantities  non-dimensionalized  using  radius  a  •  1  and  friction 


velocity  v*  «  1. 

r  (2E)7/6 

(no')1/3 

.c1 

• 

“h 

Table  C-l 

Calc.  • 

Eq.  (1.3-18) 

Eq.  (2.8-14) 

Mens. 

Est 

0 

1.78 

2.94 

3.70 

1.43 

1.0 

2.0 

0.1 

2.00 

2.98 

3.70 

1.61 

0.2 

2.45 

3.09 

4.58 

1.65 

1.2 

2.5 

0.3 

2.95 

3.29 

6.71 

1.45 

0.4 

3.45 

3.62 

7.35 

1.71 

2.0 

4.0 

0.5 

4.30 

4.16 

7.40 

2.42 

0.6 

5.30 

5.03 

7.40 

3.61 

3.0 

6.0 

0.7 

6.70 

6.52 

7.40 

5.91 

0.8 

8.20 

9.58 

7.40 

10.60 

7.0 

13.5 

0.9 

10.00 

18.80 

7.40 

25.40 

15. C 

25.0 

Notes: 

1.  Laufer 

regards  his  measured  dissipation  data 

as  too  low 

and 

attempts  to  estimate  better  balues  as  indicated. 


1/3 

2.  The  quantity  (00')  vas  obtained  from  the  theoretical 

velocity  profile  based  on  the  Prandtl-Nikuradse  mixing  length 
theory  using  Eq.  (C-l)  Refer  also  to  Eq.  (1.3-7). 


Ill 


o 


o 


I  .  a  Laufer's  Measured  Values 

0  Laufer's  Adjusted  Estimates 
Heuristic  Model 

^  ,.i  .1  . >  ■  . .  .  . . . . .  — .1.— 

0  0.2  04  06  0.8  1.0 

Relative  Radius  1 
R 


RGURE  C-l  .  INVERSE  DISSIPATION  FUNCTION 

IN  A  PIPE 


JL 

ft 
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TABLE  C-3 


AXIAL  VELOCITIES  IN  A  COLD  TURBULENT  JET 
AT  AXIAL  STATION  |  -  20 

Measurement  by  Corrsin  and  Uberoi,  Ref.  (7) 


V  -  mean  axial  velocity 

Vq  ■  axial  velocity  at  axis 

v  -  root  mean  square  axial  perturbation 

V  1 

b ^ ■  radius  to  point  where  y—  ■  y 

z  ■  axial  coordinate 

D  *  diameter  of  jet 


0 

1.0 

0.205 

0.2 

0.95 

0.21 

0.4 

0.85 

0.215 

0.6 

0.74 

0.215 

0.8 

0.625 

0.20 

1.0 

0.500 

0.18 

1.2 

0.40 

0.165 

1.4 

0.32 

0.145 

1.6 

0.25 

0.12 

1.8 

0.18 

0.09 

2.0 

0.135 

0.07 

2.2 

0.08 

0.045 

2.4 

0.055 

0.03 

2.6 

0.04 

0.02 
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TABLE  C-4 


AXIAL  AND  RADIAL  VELOCITIES  IN  A  HEATED  TURBULENT  JET 
AT  AXIAL  STATION  |  -  15 

Measurements  by  Corrsln  and  Uberol,  Ref.  (7) 


u  ■  root  mean  square  radial  perturbation 
Other  notation  as  In  Table  C-3 


0 

1.0 

0.17 

0.175 

0.2 

0.97 

0.17 

0.17 

0.4 

0.88 

0.17 

0.16 

0.6 

0.76 

0.17 

0.14 

0.8 

0.62 

0.17 

0.12 

1.0 

0.50 

0.15 

0.09 

1.2 

0.38 

0.13 

0.075 

1.4 

0.28 

0.10 

0.05 

1.6 

0,205 

0.075 

0.035 

1.8 

0.145 

0.05 

0.022 

2.0 

0.10 

0.03 

0.013 

2.2 

0.055 

0.02 

0.005 

2.4 

0.025 

0.01 

0 
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TABLE  C-5 


AXIAL  VELOCITIES  IN  A  TURBULENT  JET 
AT  TWO  AXIAL  STATIOMS 


Measurements  by  Laurence,  Ref.  (8) 
Notation  as  in  Tables  C-3  and  C-4 


■  n 


0 


0.5 

1.0 

1.5 

2.0 

2.5 
3.0 


3.5 


4.0 


4.5 


5.0 


6.0 


7.0 


8.0 


[fl 

V 

% 

V 

fV 

IV 
v  oJ 

< 

o 

V 

oJ 

k) 

0.43 

0.40 

0.31 

0.22 

0.15 

0.07 


0.082 

0.082 

0.082 

0.084 

0.089 

0.091 

0.090 

0.084 

0.074 

0.061 

0.047 


0.36 

0.34 

0.30 

0.25 

0.19 

0.133 

0.08 

0.025 

0 


0.073 

0.073 

0.073 

0.074 

0.075 

0.075 

0.075 

0.075 

0.075 

O.C72 

0.067 
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